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1. INTRODUCTION 


This report deals with the integration of two types of physiological system 
simulations. These types are classified as long-term and short-term. The long- 
term model is a circulatory system model which simulates, long-term blood flow 
variations and compartmental fluid shifts. . (1) The short-term models simulate 
transient phenomena of the respiratory, thermoregulatory, and pulsatile 
cardiovascular systems as they respond to stimuli such as LBNP, exercise, and 
environmental gaseous variations. (2-4) An overview of the interfacing approach 
is described in Section 2. Detailed descriptions of the variable interface for 
•long-term to short-term and between- the three short-term models are given in 
succeeding sections of this report. 

In order to fulfill the objectives of the study each system was carefully 
analyzed. Types of inputs and simulation forcing functions were evaluated. 

When an identical physiological variable was calculated by more than one model, 
the calculation which was most physiologically based was retained as an inter- 
facing variable. 

The major emphasis of this component of the study concentrated on the 
respiratory-pulsatile cardiovascular system with exercise playing the role of a 

i 

major stimulus. Studies of simulations involving this integrated system and 
its response to altered environmental gaseous concentrations ( 02 > 
being conducted. 
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2, OV E RALL INTEGRATED SIMULATION 
For implementatiofi of the simulation of an experiment which might 
encompass hours, days, and even weeks it is mandatory that two inter- 
facing segments be considered. One of these would handle the transfer of 
variable and parameter values during short-term simulations when all 
short-term transient models are functioning. The other interface would 
allow transfer of information in an initialization or reinitialization 
mode, .These two interfacing segments are Illustrated in Figure 1. The 
interface between the short.-term transient models and the long-term model 
is also utilized as the input for the experimental protocol. 
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Figure 1, Overall interfacing schemes of transient and 
long-term physiological system models. 
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3. IHITIALIZATIOH AND REINITIALIZATION INTERFACE 
In the following Sections brief descriptions of the salient interfacing 
features are presented. Since the interface is not individualized with 
regard to specific system types, much of the initialization of variables 
is common to all the systems. 

3.1 Circulatory to Cardiovascular 

Establishment of blood volumes after long-term simulations is 
important. In particular, the unstressed volume (V^) reflecting shifts 
due to antonomic stimulation needs to be transferred along with total blood 
volume. Changes in the resistance flow segments, such as in the renal 
component, play important roles in transient simulations of exercise. 

Whether implemented as alterations in peripheral resistance or by some 
other mechanism a cardiac output influence is necessary.. 

3.2 Circulatory to Thermoregulatory 

Skin blood flow is a variable that contributes to the particular 

compartmentalization of the thermoregulatory system, thus it is necessary 

that long-term shifts in this variable be realized. No other variables 

/ 

are strictly inherent between these two systems. That is, other 
thermoregulatory system variables are available for reinitialization through 
the cardiovascular and respiratory system components. 

3.3 Circulatory to Respiratory 

A cardiac output update is obtained from the circulatory system via 
the cardiovascular system. This transfer of variable seems logical since 
the cardiac output component is removed from the respiratory system and 
total blood flow is generated for all three short-term models by the 
cardiovascular system. 
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Long-term changes in metabolic rates must be transferred to the 
respiratory system. In addition, the blood hemoglobin (Hb) level 
variations are necessary for establishing arterial hemoglobin concentra- 
tions (Ca, , _ v) in the respiratory system. 

Cnb02.) 

3.4 Short-term to Long-term Transfer of Information 

Necessary initialization data from the thermoregulatory system include 
sTcin blood flow and a water loss variable. Skin blood flow reflects the 

Jt 

short-term thermal environmental changes as well as related physiological 
changes. Since the long-term circulatory system model does not formulate 
evaporative loss, the evaporative water loss from the thermoregulatory 
system would be utilized as an increased water loss. Consequently, the 
circulatory system would not further distinguish the total water loss and 
the evaporative loss. Details of these variable flows will be pursued in 
a future study. 

Presently, the significant variable transfer from the respiratory to 

circulatory system is the variation in Ca n cardiac output 

• (HDU2; 

influence is the dominant contributor from the cardiovascular to circulatory 
system. As the integrated system is further developed and refined, it is 
likely that other variables will be added to this initialization and 
reinitialization component. 
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4. SHORT-TERM MODEL INTERFACES 

In the top portion of Figure 1 the short-term model Interfaces are 
displayed. A closer look at the particular variables involved is given 
here. Greater emphasis is placed upon the respiratory— cardiovascular system 
interface since the immediate study concentrates on this phase. 

4.1 Respiratory-Thermoregulatory System Interface 

Only the variable which describes respiratory minute volume is directly 
transferred to the thermoregulatory system. It is an input used in 
describing water loss and heat loss formulation. Other variables of the 
thermoregulatory system which are influenced by the respiratory system are 
transferred by the cardiovascular system. These include cerebral blood 
flow and metabolic rates. No variables are passed directly from the 
therraoregulatroy to the respiratory system. 

4.2 Cardiovascular-Thermoregulatory System Interface 

There are several variables passed from the cardiovascular to the 
thermoregulatory system. Blood flows including total cardiac output, 
muscle blood flow due to exercise, and cerebral blood flow are passed to 
the therraoregulatroy system. It should be noted that the cerebral blood 
flow formulation originates in the respiratory system. In a similar 
manner metabolic rates are transferred to the thermoregulatory system via 
the cardiovascular system. Body attitude (standing, sitting, prone) as 
it relates to shunted blood flow due to physiological stress and peripheral 
resistance is transferred in a manner useful to the thermoregulatory system. 

The reverse transfer of information yields skin blood flow, a cardiac 
output influence, and a blood shunting influence due to thermal environmental 
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contributions. These are then used to update or augment existing formulations 
in the cardiovascular system. 

4.3 Respiratory-Cardiovascular System Interface 

Cerebral blood flow, described as a function of arterial CO^ and 0^ gas 
tensions in the respiratory system is passed to the cardiovascular system. 

The variable^ respiratory frequency^ is transferred to the cardiovascular 
system. Refer to Section 5 for a description of the modified version 

A 

of this expression. Instead of passing an a-v 0 ^ difference term and 
having total oxygen uptake calculated in the cardiovascular system, the 
entire development of oxygen demand is retained in the respiratory system. 
Oxygen demand is then passed to the cardiovascular system. Although not 
completely developed, arterial CO^ and 0^ tensions are passed to the 
cardiovascular system with the idea that they will be utilized in an 
implementation of CO^ and 0^ forcing for a cardiac output formulation. 

In order to fulfill the demands of the forementioned mechanism, the 
resting 0^ requirement (V02RDT) and total metabolic rate (V02DT) for a 
given exercise level is transferred to the respiratory system. The 
cardiac output subroutine is deleted from the respiratory system with 
cardiac output requirements fulfilled by a transfer from the cardiovascular 
system. 

Interface modifications are established in the following manner. The 
block diagram for controlling metabolic rate which existed in the cardio- 
vascular system is modified to the one which appears in Figure 2. A common 
interface has been established as 
COMMON/ RINTR/ROUT (10) , CIN (10) . 



Figure 2. Revised block diagram for controlling metabolic rate in the 

cardiovascular system. Variable ROUT(l) is obtained from the 
respiratory system. 
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In the respiratory model. Subroutine BC12, the output variables are 

indicated as 

ROUT(l) = AV02DM/1000 
ROUT (2) = FREQ 
ROUT (3) = C(ll) 

ROUT (4) = F(7) 

ROUT (5) = F(l) (4.1) 

where 

AV02DM = a-v O 2 dif f erence,ml02/min, 

FREQ = respiratory frequency, bpm, 

C(ll) = cerebral blood flow, 1/min, 

F(7) = PaCCO^), mmHg, and 

F(l) = PaCO^f, mmHg. 

Likewise, the input variables are indicated as 

C(10) = CIN(l) 

V02DT = GIN (2) 

V02RDT = CINO) (4.2) 

where 

C(10) = cardiac output, 1/min, 

V02DT = oxygen required for particular work load, 102 /min, and 
V02RDT = oxygen required for resting state, 102 /min. 

In the SS02W(X) Subroutine of the respiratory system, the following 

modification has been performed. The former subroutine statements were 

replaced by 

COMMON/RINTR/ROUT (10) , CIN (10) 

V02RDT = CIN (3) (4.3) 

SS02W(X) -• V02RDT - .0500 + (.0004850815 * 6.12 * X)/.25 (4.4) 

where X = work load in watts and the other terms are as previously defined. 

In Subroutine RC12 the following statements were added such that they 

appear in both the increasing and decreasing work load' paths. 


IF (WORK.LE.0.0 .AND. NWREST.LT.l) RMT(2) = CIN(3)-C(26) 
AV02DM = (F(9)*C(10)-F.(13)*(C(10)-C(11))-F(12)*C(11))*1000. 
AV02DF = AV02DM/C(10) 

ROUT(l) = AV02DM/1000. 

IF(WORK.GT.O.O) ROUT(l) « RMT(2)+C(26) 
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ROUT (2) = FREQ 
ROUT (3) = C(ll) 

ROUT (4) = F(7) 

ROUT (5) = F(l) 

The transfer of cardiac output from the cardiovascular to the 
respiratory system was handled in the main program, GRODIN, C(10) = 
CIN(l). Additions to COMMON/R/ in RC 12 include RMTM and TCT. 

Refer to Appendix 7. Lfbr the program listing illustrating the 
implementation of these statements in Section 4.3 and the corresponding 
changes in the cardiovascular system. 
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5. MODIFICATION OF I>TDIVIDUAL RESPIRATORY SYSTEM 

Improvement in the individual respiratory model was suggested in a previous 

research report. (5) These modifications were made and presently, exist in 

the latest version of the respiratory system model, (6) In addition, a-v 0^ 

difference (AV02DF) and dead space volume (DSVOL) have been added to the out- 

/ 

put routines. 

These modifications are summarized here. In the original model 
respiratory frequency (FREQ) was given by 

FREQ = 8.1 + 7.815 * (EIIT(2) + C(26)) (5,1) 

with 


RMT(2) = O2 metabolic rate of tissue and 
C(26) = ©2 metabolic rate of brain. 

Thus, Equation 5.1 didn't respond to any forcing other than 0^ demand, 
formulation was replaced by 


FREQ 


ffl + 32 / 

i + aUc - l\ 

\k \ 

a /DSVOL/ / 



(5.2) 


This 


with 


RC = 0.015 mlUj 

o 

V = expired ventilation, 

A E ' 

a = 1 95 = inspiratory elastance 
expiratory elastance 

DSVOL = dead space volume, (5) 


h. 

h 


and 


Upon substituting the constants, FREQ is given by 
FREQ = ((1. + (.726 * VE) /DSVOL)**. 5 - l.)/.363 

with 

DSVOL « 0.140 + 0.002 * VE. 


(5.3) 


(5.4) 


Dead space ventilation, originally defined as 
DEADVT = .1107 * FREQ + .0785 * VE 


(5.5) 
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is now given by 

DEADVT = 1 + .098 * VE. (3.6) 

The representation for minute volume (TVNT) remains- unchanged. 

The a-v O 2 difference expression that is necessary for the integrated 
system to function is added to RC12. Here, 

AV02DF = AV02DM/C(10) (5.7) 

where 

C(10) = cardiac output and 

AV02DM = (F(9) * C(10) - F(13) * (C(10) - C(ll)) - F(12) * C(ll))> * 1000 
as defined in Section 4. 

The terms in AV02DM are defined in exactly the same way as in the 
original respiratory system model. See Appendix D of the listed reference. ( 7 ) 
Also, in comparison of the original respiratory program and the latest 
modified version there are additions to the output statements of RC12. (5-7) 
Statement #*s 218, 219, and 220 have been modified to include AV02DF and 
DSVOL. In a similar manner statement //'s 246, 263, 264, and 265 now 
include these new output variables. 

Before leaving the discussion of the calculation of a-v O 2 difference 
the significance of transport delays should be considered. Slight errors 
exist in the present formulation. Blood flow transport delay times are 
not considered in the calculation of venous blood concentrations. See 
Equation 5.7. Actually the concentrations and compartmental blood flows 
should correspond to the same time. This means that an additional book- 
keeping operation should be implemented such that the concentration at 
the lung entrance reflects the delay times and their corresponding 


contributions . 
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6. EVALUATION OF INTEGRATED RESPIRATORY-CARDIOVASCULAR SYSTEM 
The evaluation of the Integrated respiratory-cardiovascular system 
proved quite encouraging. Several types of simulations were tried. A 
200-watt exercise level of 5 minute duration was used as the stimulus. In 
this section two systems are compared. The basic difference between the 
two systems involves the- formulation of metabolic requirements. System A 
is shown in Figure 3 and System B is shown in Figure 4. Selected responses 
for these two systems are illustrated in Figures 5-14, 

The major variations in the responses can be summarized as follows. 
System B doesn't allow for the rapid increase in heart rate that occurs 
with System A. Since the cardiac output doesn't vary appreciably between 
the two systems, the over response in heart rate of System A is accompanied 
by a decrease in stroke volume. System B is a slightly more efficient 
system since a lesser amount of 0^ (0.1 102 /mln) is required to sustain the 
simulation at this steady-state exercise level. The differences in the pre- 
exercise variable levels are related to the differences of basal conditions 
for the respiratory system and resting conditions for the cardiovascular 
system. This feature is coupled with the fact that the- cardiac output 
settles to5V 6.8 1/min compared to the 6 1/min for the original respiratory 
system. After considering all of the variables' responses and the control 
of regulation involved with each one, System B seems to perform in a more 
satisfying manner. Therefore, the system shown in Figure 4 is recommended 
as the Integrated system for exercise simulations. 









15 


Pulsatile 

Cardiovascular 

System 

Cerebral Blood Flow 

Transient 

Respiratory 

System 


Respiratory Frequency 

■ 


Transient 0^ Uptake 



Cardiac Output 



Resting 0., Requirement 


i| 

Steady-state 0^ Requirement 






Figure 4. Respiratory-cardiovascular system Interface which utilizes the 
metabolic formulation of the respiratory system during exercise 
stimulation. 
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01 2 345 67-89 10 11 

time ,min. 

Figure 5. Inspired ventilation versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response. (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements, (b) Metabolic requirements controlled by respiratory system model. 



Respiratory frequency, .bpm 



Figure 6. Respiratory frequency versus time for five minutes of 200 watt exercise stimultion and corresponding 
off-transient response, (a) Individual cardiovascular and respiratory system simulations of 


metabolic requirements. (b) Metabolic requirements controlled by respiratory system model. 



Heart rate, bpm 


210 



time, min. m 

Figure 7. Heart rate versus time for five minutes of 200 watt exercise stimulation and corresponding off- 
transient response. (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements, (b) Metabolic requirements controlled by respiratory system model. 




Cerebral blood flow, 1/min 



time, min. 

Figure 8. Cerebral blood flow versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response- (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements. (b) Metabolic requirements controlled by respiratory system model. 
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32 
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Figure 9. Cardiac output versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response. (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements. (b) Metabolic requirements controlled by respiratory system model 




stroke volume 
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Figure 10. Stroke volume versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response, (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements. (b) Metabolic requirements controlled by respiratory system model. 



Tissue 0 tension, mmHg 



time, min. _ 

Figure 11 • Tissue O2 tension versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response. (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements. (b) Metabolic requirements controlled by respiratory system model 
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Figure 12. Arterial 0^ tension versus time for five minutes of 200 watt exercise stimulation and corresponding 
off-transient response. (a) Individual cardiovascular and respiratory system simulations of 
metabolic requirements. (b) Metabolic requirements controlled by respiratory system model. 
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40 
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time, min. 

Figure 13. Arterial-venous 0^ difference versus time for five minutes of 200 watt exercise stimulation and 

corresponding off-transient response, (a) Individual cardiovascular and respiratory system Simula 
tions of metabolic requirements, (b) Metabolic requirements controlled by respiratory system model 
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time, min. 

Figure 14. Tissue 0^ metabolic rate versus time for five minutes of 200 watt exercise stimulation and 

corresponding off-transient response, (a) Individual cardiovascular and respiratory system simula- 
tions of metabolic requirements. (b) Metabolic requirements controlled by respiratory system model. 
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7, APPENDIX 


7.1 ProRram Listing of Particular Subroutines 

Computer 'Program Listing of the subroutines of the respiratory system and 
pulsatile cardiovascular system models which were modified for the integration 
of the" two systems. 
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99 C cerebral blood PLO.V. 

lOQ C - 11 03 - - — 

101 C TISSUE BLOOO FLOW. 

- 132 C 12 OT 

103 C ARTERIAL H. CONCENTRATION. 

108 C I 3 CA<H*» 

105 C ARTERIAL 02 TENSION. 

106 . - - C - 18 -- PA(02) • - 

107 C 

lOa C IG 

109 C TOTAL gas C ONCENTR A T I O.jS AT BRAIN EXIT. 

no c_ . 16 . cvB(Co 2 i cvaioji *• cvb(n 2 > 

ni C total GAS CONCENTRAT IO«S AT TISSUE EXIT. 

112 C- 17 — CVT(Co2»-»— CVTiO?) * CVT«.\z) 

113 C TIME. 





T 


3B 

c 

13 

37- 

c — 

ARTERIAL 

MO 

c 

14 

H 1 

‘ - c 

arterial 

42 

c 

IS 

43 

- -- c- 

total Ga 

44 

c 

16 


C 

C 

c 

- - c • 

c 

c — 

c 

c 

c 

c - 
c 

c 


VTRANUsI 

arterial gas concentraTioms aT brain entrance. 

I CA8<C02) ■ CA(C0?)(T - TAB) / 

2- -CaB(02j . CA(02),T - TAB) 

3 “ ca‘n2>‘T ! 

VENOUS BRAIN 6A5 CONCENTRATION AT lUNG ENTRANCE. 

H • CV3<C02KT - TV8) - - 

S CVB(o 2> *T - TVd) 

A CV 81 N 2 ) IT - TVa) • 

— VCNOOS— tissue GAS concentration -aT' LUNG--ENTRaNCEi 

7 CVTtCo2) (T - TVTl 

— - 8 CVTI02) IT - TVT) ’ - ' ' 

9 CVT«N2) (T - TVT> 

ARTERIAL GAS concentrations AT TISSUE ENTRANCE. 

10-- CATIC02I - CA<C02)(T - TaT» 

11 CATtOEl ■ CA1021IT - TATI 

--12 CATIN2) ■ CAIN2MT - TATI - ■ - 

-ARTERIAL rt* CONCENTRATION AT CARQTIO 300 1 E S' 5 I TE — = 

13 CAO(H>) ■ CA(H+)(T - TAO) 

"ARTERIAL 02 TENSION aT CaROTIO BOO I £ S * S I TE . T" 

IN PA0I021 ■ PAt02)tT - TAO) 

arterial h+ Concentration at brain entrance. - -■ - 

IS CA3(N-.I ■ CA|H+)(T . TAB) 

total Gas concentration frou bRain at lung' entrance. • ■ 

ti (CV3IC02) ♦ CVB(n2» ♦ CVS(N2))tT - TVS) 

-Total-saS concentration from tissue at Lung entrance.-- 

17 <CVTCC02) ♦ CVT<n2) + CVTIn2)HT - TVT) 

0( IBl 

• FOR ons) THE SYMBOLS BbBAROmETRIC PRESSURE, < 17 »WATER VaPOR PRESS.^- 
K«CONVERSION Factor for .atm to MMHG. A«S0LU3ILITY coeff.of gases, 

— HbCOMPOTER time-step-, HBbBLOOO- oxygen -capaci ty 




o s 

■■©'Sh 

?o r- 


C 3» 


AN2 <3 M7> 

A02 13 — 47) 

An2 (B - 47) • 

>16 * 2.3 IH3 ) 


60 
61 
62 
A ^ — 

C 

• -c 
c 

■ 7 

- • - 10 
1 1 
1 ? 

863/(3 - 47) 

0.62 • — - - - 

K ACSF(C02) 



64 

c 

is 

K aCSF(N2) 


6S 

_c 

14 

2»H ■ ■ 

. . 

66 

c 

IS 

1 .97»H 


67 

c 

FI20J 



63 

c 

Comp artmental gas tensions ano Concentrations. 


69 -- 

- “C 

1 - 

PA(02) 


70 

c 

2 

< AC02 Pa(C02> 



LA 

NA 



17l 

C 

3 

P6 102) 

t ?2"*" 


H • 

*: AC02 P6IC02J 

173 

c 

5 

PTU)2> 

1 7H-— 

— c 

6 ■ 

K AC02 PT(C02) 

175 

c 

7 

PA(C02) 

17* 

c 

• 8 

PAI021 

177 

c 

9 

CA(q2) 

178 - 

— c 

— 10- 

CAIN2) • 

179- 

c 

1 1 

CA(C02) + CA102) * CAIN2) 

—-180 - 

- c 

■ 12 

cva (02 > - - ■ 

181 

c 

13 

CVT(02) 

162 

c 



183 

c 

product 

OF DIFFUSION COEfFS.ANO GAS DIFFERENTIALS ACROSS dLOOO-BRAlN 

1 8H 

c- 

—barrier 

- - ... 

135 

c 

IH 

DC02 (PB(C02> - PCSF(C02)) 

186 

- c 

IS 

002 (PB102) - PCSF102)) 

187 

c 

1* 

0N2 tPBtN2) - PCsF(N2)) 

- - 1 63 

• c 

. — 

t « 

189 

c 


P3(02> 

l?o - 

C- 

— • i'a 

P9(n2> - • ... - 

191 


OIMENSIO.N XNaiH.2). O.lIH), lDjl2) 

192 


COMHON/2/ C* XH. SV, vTRAH, RK( SCt OCt Af 0» F, VOUt RMT> BC* <JF» 

193 


1 

TaU» CC, CnB, CH, CPH . DO. VE. VI, CP8 , CPT. CADK, X, OT, 

' • 19H 

- — A . 

2 

-.IRK, LOC. ITERx, INDEX, ]. J, rt, N 

195 


COMHON/R/ XoSiXhH.CXT.WoRK.DOmI ,DU rt2>DUn3iY(0RK2iKMTaf R mTb2*TIh£0F 

- • -1 96 • 

— 

1 .SHLlN.lTTY - 

197 

c 

ITTY ■flg for tty MOoE* 

193 • 

c- 

— o« output To prIntEr oatch Hooe)« 

l99 

c 

• TTY 

TTV I/O AND 1st TI.HE TO SUBROUTINE RClZ. 

.. 200 • 

--C 

1 > 

TTY I/O ANO NOT I ST TtHE TO kC12. 

231 


DATA 

ITTTY/*TTV •/ 

202 — 

c* 

data 

FOR Initial conoitiohs • -- - 

233 


C( 10)«CtNU) 

20H- 

— 

-IF (CXT.GT.O.) GO TO AO 


205 . flSITE Ui3> 

206 S FORrtAT t/* - GROOINSt HESPIRATOhV CONTROL ■ KODEl. V/ ) 

207 300 CONTIUUe 

'2oa--- »-tRtTEJ4.R83> ... 

209 S83 FORrtATMOADO DATA...’) 

210 -- C-- reao indication oc Batch or ttt mooe» 

211 REAOtS.HSO) ITTT 


212 

213 

El-F- 

215 

- 216 
217 

213 

.J 220- 

22 1 
222 
223 
228 
225 

22* 

227 


880 FoRrtAT ( AH) 

IFIITTY ,NE. ITTTYl ItTY-« 0 
rtRlT£<4,90) 

90 FORHAT I Irtl ,IX.37H»RE<;P1RAT0RY CHEMOSTaT — INPUT DATA*/) 
-■c ■ data for initial CONOtTIONS I 

00 i'O I B UNO 

C IlOb HAS PPOaLE.’l ..ITrt ENO** « SO TH J S ISUT USED TO 
C OtTER.'UNE end Of HOHlNOn capability To 5T«KT another 
-C- • hOOEL Ron in same computer run). 

REA015>I?0<EnO~30|) 

10 CONTINUE 

C establish COntJUTER STEP IHOEPENoEnT OF JNPUT (JATAt 
C(a*)«.78125E-2 
l9o FORMAT I 5 A , F I S . 0 .5 A . 2 A H > 

..DO 20 I * 1 *H 

IPsO ■ I * HO 



<xz^ 

22s^Tg>- 

~5 ^ 

■< -t'? 



228 

read (5>190) «C(I)t 

(XN8( I »J) 

“229 20 

continue ■ 


230 

00 30 1 ■ 1.2 


- 231 

KEAO (5.I90> rntu.)* 

(XNM I .J) 

232 

IP10 » 1 ♦ 11 


- 233 30 

CONTiHUe ■ - 


231 

00 10 I ■ 1 .2 


- 235 

READ t5.190)-0J(I>* 

-(XN8( I » J) 

236 

IP1Q • 1 ♦ 16 


237 -- - 10 

continue - ' ' 


238 C 

239 C OUTPUT INPuT DATA. 


210 

J * 1 



J B U2> 

J “ 1»2) 

J “ - •- 


2‘Jf — DO 75 I -- 

2HZ JX > J H 

213 — - aKlT|-(5.?2) J,(cnZ)f T2«J«JX> 

211 92 FORHATC* » . I2.2X.5tF9.1) ) 

~ ■ 21 5 J ■ J ♦ 5 - * 

216 75 CONTItlUC 

217 — t(niTE(6,92> -Jr<BC( I > H-J .1) - 

216 J s IS 

-- 219 - — n»IT£(6.92) J,RMTU>*RMT« 2 > tOJt J> ,DJ(2> 

250 C 

2Sl — C IF TTY I/O HAX.TIME rtll-U COME FROM WORK CARD, -- 

252 IF(ITTV .NE. 0> C(IS) >« 9999999999. 

251 C FUCOZ) 

" ' 255 OUH I^C ( 31 > 

256 C Fl(aZ) 

- • 257 0UH2»C<321 

258 ■ C FI 1142) 

259 0U«3BC{33> — - • 

260 *OHX“0. 

261 -- - - H0R<'2b0. 

262 C META30LIC RATE OF 02 CONSUMPTION IN TISSUE. 

“■ 263 R«TB“C1nI3) -C (26) 

261 «MTB2>»CIN(3)-C{26) 

266 TIME0r»0. 

- • 267 -- ••• X01 = 0. ■ - •• • ■ • . 

268 XMHBl0.»C(36)/0.0078l75 

269 ■ ' MMI1“6 

270 201 CONTINUE 

. 271 XOS«XOS-»XMH — • - 

272 IFIMMM.CO. 1 )XOSbXOS6C(36) 

273 HMri-1 . 

271 C(35)«<0* 

275 C(10)«0. 

276 C 

™ 277-- C INITJAL GUESSES FOR ITERATIVE LOOPS 

278 C ARTERIAL CO.NCENTRaT I ON OP C02. 

... 279 - — -• — CC« 1 > ■ 0.6 

280 C BRAIN COnCENTRaT-ION OF Co2. 

281 ^ CCC2) - C(1) 

232 c TISSUE CO'^CENTRATION OF C02. 

2S3 --CCiS) • C(7) 

281 C BRAIN C02 tension. 






(SS)"!® 



— - 





28S CpB ■ SfJtO 

2nb c --T issue coz te«sion»- 

237 • CRT ■ 50*0 

238- lF{JtDS.GT«XHH) G0T0207 

239 c SETS VASJOUS CONSTANTS AND aGGReGaTeS OF CONSTANTS 

290- - ■ - C- TmAX. 

291 C( IS) * Ct ISl ♦ .0001 

292 C-- PRINT AUL-TIHE, — 

293 C C(39) ■ C(39) ♦ .0001 

— 299- --C Factor of i-e-7 HUtTiPuYiNG oiRFiisjoi.- coefficients. 

295 00 200 I “ 27»29 

296 C(l)»CtU»l .e-7 

297 200 CONTINUE 

298 - - -202 — CONT INUE ■ 

299 IRX « I 

- - 300 - — - 

30 I N a 5 

30 2 lOdt I ) a 0 - - — 

303 C solubility COEFFICIENTS. 

304 C All )a. {alpha )C02.- A < 2»a-^ alpha ) 02 . At3)« <ALPHA)N2» 

30S C A(9)« <ALPHA)C02i A(S)a (AlPHA) 02, A(6)a (ALPHA)NZ 

306 “ O.Sl _ , 

307 ^ AI2) a O.024 

- 308 -- -- A13) a 0.013 -- . 

309 AIR) a 0.51 

310 AlSl a D.02M • - •- - - 

311 A(6) a ,0»013 

" ■ 312 C- ATrt/HrtHG CONVERSION FACTnR. 

313 SK « 0.00132 

3ji< CASdONlC acid dissociation CONSTANT. 

3tS CADK a 795.0 ' • 

316 C VOL< I >-VOL ( 10 >a VOtUriES-USEO IN CALCULATION OF -V aR 1 ABLE- T 1 «E- DELAYS.- 

317 VOL ( U = 0.015 

- 313 - Y0L<2) a 1 .092- — - 

3 1 9 VOL (3) a 0. I 83 

-- 320 VOLtR) a O.Oo • ■ - - • ■ 

321 VOLtS) a 0.183 

-32 2 VOL 1.6) a 2.9M — . . _ 

323 VOLt?) a 0.73s 

--- -329 VQLlai a 1 .062 

325 V0UI9) ■ O.QOa 

-326 V0L-(10)a l,0o2- 

327 C 

328- c (METABOLIC RATe OF C02 IN_8RAIN TISSUE. I / SAH£ FOR 02 

329 qF( 6) a (C( 2 S> •* RMTtl) )/(C(26) ♦ RHT(2)1 

330— --C — B-47 • ™" — - - 

331 0( 1)»C«3C)-47. 

— - 332 - - 00 2 10 I " 2.4 

333 C products OF COfIVEHSION FACTORS AND SOLUBILITY COEFF I C 1 Ei'.TS » 

-■ 334 0(t) a 5<«A<1-U - • 

335 0( 1»9) a SK*A( 1*2) 

336 C — 

337 0( 1*3) a 0{ 1 1 aon ) 

33ft 210 CONTINUE 

339 C FACTOR USED IN ESTABLISHING ca(C02) 

— 340 — - ■ ■■ 0(31 a- 0.16- a.— 2.3.C1 17) 

34 1 



C 



3*»2 Of * ■ ai3.0>'0 ! I 1 

3H3 c — F actor used in ESTA»LiSHtN<; cstco2>* 

3H4 D(iO) ■ 0.62 

3^s — t -Manipulation op {^ohpuTeR Ti«e st^f. 

3<)6 D< t'l) ■ C(36)«2»0 

3<(7 * -0(j5> ■ D<iH> ~ ♦0l*C<36) 

3se c 

3*19 CALL- Rc3 

350 call Rc9 

--•351- CALL RC5 tCp8. F < <♦ ) , C ( <♦ J , eC<2J) 

'352 CALL RC21 {CHB{2>* F(3l« FIR}. C t H ) i CHU). CPH( 

•• 353 • • call RCi9 <CPa* CHBC2). CC12>* ftC t 1 ) » F I <* ) ) 

35** CALL RCS (CPT, F»6), r(7>, BC(3>) 

35S- CALL RC2 1-- <CHB< 3> » F ( *4 > . F(6». C<7). CH<3>i CPH{ 

356 CALL HC19. (CPT, CHt* 1 3 ) , CC ( 3 ) , gC ( 1 } • F(6)) 

-• 3S7 CALL P.CPO - -- 

35a CALL KC7 

359 CALL RCa -- — ■ • 

360 CALL RC9 

361 call-rcio — - 

3b2 call RC 1 1 

— 363 - ■ • CALL RC 1 2 - • 

36*4 GO TO 60 

- • 365 -- SO CAlL RCiS 

366 CALL HCJ6 

367 =~60 Call RCI3-— ■ 

368 CALL RC12 

■ 369 - • -C - . ... 

370 IF1C(35).G£.XMHJ 60 To 20l 

T 371 " C - 

372 IP <C<35) .GT. C*ISH G0T080 

373 — ■ - 1FICXT«6T»C( IS) » -GOTO 80 _ 

37<* 70 call RCIH 

'-■•375 •■-UU»AK00lCt35)»Dtl‘U) 

376 IF (UU .LT> ,0001 ,0R. UU ,CT« 0(15)} SOTOSO 

-- 377 - - - “ — RETURN .... ... 

378 C 60 TO 60 

•—-379 So «RtTt»6.7S» - .... 

330 78 FORlSATt*! FINAL VALUES FOR FOLLOWING VARIABLES*' 

331 - •- .IF (C(37) ,6T. l»0E-5) GO To 2&0 

352 220 CT£R)4 a 0.0 

383- IF (VTSAN(l*r} - lOA.O) 230, 2*(0, 240 -• ■ 

3dH 230 CTERN « ( 2 3 * 6E-9 ) • ( ( 1 o') « O - V Tft AN ( I *)))••**• 9 ) 

385 2So-C(37> « C ( 20 ) *• I C U 6 ) • VTR AN ( 1 5 ) + {1.0 « C<16))»C 

336 I + C ( 2 1 ) *VTRAM ( I 3) ♦ CTERN - VI 

387 I ■ 37 ... 

333 aRITE ( 6. 192) I ,C< i > , (XNU.J), U > U2) 

389 250 00 260 I » 1 , IH- 

390 aS!TE(6,IV2) 1 ,CII ) , (xNJt.J), J •• 1,2) 

— 36J 260 continue - 

392 xfilTE {6,19*4) 

— 393 6R1T£{6,S30) 

39*) 630 FORMAT t '0)*ORHAL TERMATIOM') 

- 395 301 CONTINUE 

396 STOP 

397 C ■90-ForHAT * lHl*(8x37H,ftE56lRAT()RY CwEnOSTAT — • iNpUT 

398 C 92 F0){MAT •{‘*2X13, 10XF10.6. 10X2A6) 




C'>90 FORMAT (BXF I S .0 .<>AZA6 ) 

192 format I •-’•.l»,2A.FlS..*;,2X.2A‘t> 

19R FORMAT (IHl) 

ENO-- - 



036-003‘»32*TPFS.RC12 



- 5- 
6 

• 7- 

a 

■ 9 

to 
1 1 
12 
13 

1 ** 
15- 
lA 
17 
le 

•i9 

20 

2 1 ■ 


— SUSROUTIME RC12 - 

CO!IMOiJ/«lHTR/ROUT( 10) ,C1N< lOI 

-r 0lH£rtS10!l Ct*iO), XNlsnt2)t S(/(l8»b0)* VTRANUftli RKIMtH* - — 

** 1 50*14.5)1 Dc'l**’. A**’)* FC2 Q)i VOLUO). RHTtZ). 

2 BC(4). tAo15)i CC<3). CHB«3>i CH(4)« CpH<3». 

3 00»4) 

COMMOM/Z/ C» XH.-SV.- vTRaN. RK. Sc . OC t A » -Ot F* VOU. RMT-. -BC « OF 

J TAU, CC. Cri 3 , CH. CRH. t)Q, VC, VI, CPB, CPT. CaOK, X, OT, 

- - 2 IRK. LOG. ITEkV, INDEX, t, J, H, N 

COMHON/H/ XOS,Xrtri.CXT.AOHK ,OU:il , DUf 12 , OU M 3 t 1 ( 0 RK 2 , BM TD , K. 4 Ta 2 » T I HE OF 

1 .RKL IN, I TTY , ITTYOT » t T TV IN (HRKTTY (50, 3 ) , 1.E XE C . H Afi KER , NWRES T . ... 

2 .RMTM.TCT 

Of,XA IRUM/ifiON »/, ISToP/' 5 TOP«/,HORe/««ORE*/- ■ • -• 

DATA iback/ •Back*/’ 

c- ■ DIMENSION rtRXTTY (50,3) 

C6969 FCKMATIIh SMSUB RC12) 

■ -C- - OUTPUT PUNCHED CAfinS AND PHInTED 

CXT»C(35>-.XDS-10. 

— IF (CXT.I.E > 0 . >CXT“+G. 

c dead space volume 

05V0U“0. J 40+O.00Z»VE 


22 C RESPIRATORY FREQUENCY. 

23 — FREJ=( ( I .* ( .726«VE ) /OsVOL) ••»5"l • > X.3A3 -- 


24 c dead space VEnT.IlaTIoN 

25 DEAOVT" I .-^.OYB'VE - ■ , . . 

26 c C < 3 1 ) = { OE ADVT »C t 1 ) ■►V£*DUM 1 ) / ( oEADVT'f VE 1 

. .. _ 27 — -C C(32)«CD£AOVT«CC2)*V£.OUM2)/(OEaOVT*VE) -■ 

2a C C(33)“<OEADVT«C(3I*VE*0UH3)/(0 EaDVT+VE) 

•- 29 c minute volume. ■ - - ' 

30 TVNT»DEA0VT+(v£-t’Vn/2. 

RATE. - . . 

32 hRAT£<.43.a«(RHT(2)*C(2in*5 4.5 

33 C • • 

34 C 

- - 35 - - IF(CXT .LT. TIHEOF) Go TO 203 • -• 


3e C 

37 C HERE- tF NEED- TO READ A NF.V WORK LOaO-CARO. • 

3B C branch if in batch MODE. 

39 - - IFUTTV .EO. 0) GO 70 SOO 

40 C 

42 C HERE IF TTY MODE.... 

N3 — — IFUTTY .EO. 1) GO TO 550 - --- • 

44 C HERE IF TTY MODE. AMD 1ST TIME THIS ROUTINE CALLED. 

MS - ITTT o 1 

4A nRITE(6.S0Sl 

47- • SOS formai ( »o I nput iVORX Cards...*/) 


43 

C 

1 

» 

rtORXn 

OHX LOAD 1 UaTTS ) . » ♦ ♦ / 

.. m9— 

-c- 

2 

1 

Ml NS* T 

IME FOR »0f?< load...*/ 

5D 

c 

3 

# 

print* 

TIME INCR1m£NT(«INS)F0R PRINTOUT..,*/ 

• SI ■ 

■ -c 

...M 

9 

EXEC • . . 

•/ 

52 

c 

5 

f 

MOH£e 

INPUT MORE before EXEC...’/ 

- 53 

-c 

- 6 


run « 

EXEC.AITm A6DVE.THEN CAN INPUT AGAIN..*'/ 

54 

c 

7 

f 

ST0P4 

EAEC.AlTh above TrtEN STOP...’/ 

-55 • 

-- -c 

a 

1 

BACK* 

erase PREVIOUS WORK rECOrO...*) 

53 


504 ITTYIN ■ 0 




S7 

sa 
s? 
-60 
61 
62 
63 
6*1- 
6& 
66 
67 
- -68 
69 
— 70 
71 
-72 
73 
7*. 

75 

76 

77 


ITTYOT n I 

c . . 

SCI IFIITTYIN .Lit SO) 60 TO SC6 

- — c Here, if buffer for *iO«r i oao cards ts fouu. 

6R I TE (6 ,S1 1 I 

- - 511 FORIIaT ( 'CftOFFlFR FOR V.oRK LOAD RECORDS FULL.'/ 

1 * 61LL USE EXEC" Run. * 1 

• LEXEC . 1 RUR - 

GO TO 5Sl 

- C 

50 6 ITTYlN n ITTYIO * 1 

.. — 509 ftRjT6(6.507> 

507 FORilATli .YORK HINS PRI.JT tXEC 'i 
1 '(F6.2,tX.F6.2.1XtF6.2.1X.AR)...') 

RE ADIS, 502, ERR” S09J (rt RKtTYUTTYIII,JJ,J« 1,3),LEXE.C 
_ 502 format <r6..? , 1X.F6.2, lx ,F6.2, I X, A*U 

6ftlTEl6,50 3) (.•.RXTTYliTTylH.J),J‘»l, 31, LEXEC 

503 Fr)RMAT(3(i i,F6.2J,A<i) 

IFlLEXec .RE. IBACKI RO TO 518 

ITTYIN- o ITTYIH - 1 

IFUTTYlfl ,lT. 1) ITTvIn m 1 


7a GO TO 509 

79 C 

»0 518 IFCLEXEC •E-'7* IRLUJ .on. LEXEC -ca. ISTOR) GO TO 55 I 

31 1f(l£XEC .£0. hOREl Go TO &01 

S2 wR 1 TE ( 6 ,5 10) - 

83 510 FORMATl' E XEC .P aR A ME TfR .(RONl, TRY ACAlfL'l 

39 Go TO 509 

35 C 

S6 - - c HERE IF 1ST TIhE THIS KOoTtr,E CALLEU, 

87 C SEE IF MORE YIORK CARDS IN BUFF ER 1 ,’IR tCT T Y 1 SO 0 , 3 ) ) 

-SS SSa IFllTTYOT .uE.-lTTYlN) go TO 55 1 . 

39 C HErtE IF EXAJSTED .'.'ORX CArD. BUFFER ( ..R K T t Y t bOO , 3 11 . 

93 . 1F(LEX£c .EG, lAUi.) Go TO 509 

91 C FORCE END Or COMPUTER RUN ’AhEN LEXEC“ 'STOP*. 

- 92 Ct 151 - 0. - 

93 GO TO 1210 



95 SSI .10RK2 « NRKTTYt ITTTOT. 1 1 

96 OURAT a *KX T T Y U T T T , 2 1 . - 

97 C(39) s ARXTTYI ITTYOT.3) 

93 . . . iT-roT s ITTyoT ♦ I 

99 GO TO 606 

100 c 

10 1 C 

102 C-. . 

103 203 IF I MARKER. EO .01 GOTOlfil 

ICR - - 1 . )i0RK“i.'0RK2 

105 .MARKLRa] 

103 .C--STSTEM KESPOUSES: Tl/iE Cnt.STA’JTS FOR f-ORK LOAD LE VtLS 1 1 NCKE A S 1 N G 1 , 

107 IF ( NCRK.LE.O. 1GOT02 

ltl8 IF < uORK .GE.SO. 1 TCT = 2. 3/ ( 2.* AORK/200. I 

109 lF(AORK.LT.SO)TCTn<,.6 

110 C tissue 02 KETA30LIC RATE. 

111 KH T|2j»'>S 02.V( A OR <1-15502. V( .0«KI-KHT321»LXP(-TCT»(CX1 -TIMEOMll 

U2 VT IH£»1 .1-1 . l«£X?l-TCT»lcxT-TIMeOU/l.V2) 

113 C term used In VI that is a CDMPOuEIjT OF transient RESPONSE RELATED 



TO WORK t-OAO* 

RmlIH «SSO2WtwORK>-»Ss02.^<v.0RK>-R(!Td21 *U4-VT1«E1 

IF t VTI kE.GE. 1 . > RHLlNaSS02W(«0KK> 

•TiSbUE C02 «ET&aOLIC RATE, 

RMT< I l».86«HrtT(2> 

IF jC (35 1 .LT .C{RO j ) G0T02 

wRlTg (6,333) RHT( } ) ,RHTi 2>- - . . 

J3 FORMAT) » 0 * , lX,2bHCrtAMGE IN METAbOuIC K ATES , SX t 7 RCoZ" 

J 5K,6rtrt«02» ,FlO,M,/| 


GO TO 


2 CONTINUE 

!E < f.ORK«l.E,0.0 •AN04 _NWReST»LT. n RHT ( 2)»C IN<3) 

AV0 20MA(Ft9>*CUO)-r(73)«(C<lO)-C(ll))-F(l2)»CU 
AV02 UF»>aV020“./C( 1C I 

R T ( 2 ) ” F K L 0 ~ I F (*AoKK . 6T. O.O) f?OUT(l) = RmT(2.)4 

ROUT (3) =C 111) 

ROUT »aF { 7 » ' -= — . . 

ROUT (5)*F ( 1 ) 

“ U = AHOO<C( 35) , 0>5) 

IF (y .t,t» l,0E-5 .OR. U ,GT* •4999) GO TO 

IFlC(35l.yT.C(40) )G0To)230, 

C(40)'-C(40)4C<37) 

■ARTERIAL N2 TENSION.-"- •• ■ - 

1210 PAN2 o (X I ) «C( 3 ) 

• TISSUE 02 tension. 

PT02 « C t 3 ) /y ( 3) 

Tissue n2 tension. 

PTN2 n C (9 I /O ( 4 ) 

-cerebrospinal fluid PH ♦ EQUATION 6»2 • — . . 

PrtCSf a 9. ~ RCFHCH14)) 

VENOUS cRaiN ri+ concentration , EQUATION 4.7 « - - 

• RV3 » CAD<*F(4)/tCC(2) - F(4)> 
venous brain PN , equation 4.6 • . ■ 

PHV0 = 9, - RCFKHVB) 

■•■VENOUS Tissue h'» concentration , equation 5.7 - 

rtVT = CAO<C‘F ( A)/ ICC{ 3 ) - F(6)l 
VENOUS TISSUE PH , ESUaTtON S.A . 

PM7T ■ 9. - PCFUHVT) 

-RESPIRATORY QUOTIENT (AwvEOUAR). 

S-3 » I (C» U > •VTRAI.,(4) ♦ OFl l)*VTRAN(7) >/CU0) - 

1-.- (F(9) (Ctll>*VT-RAN(5l-+ OF ( 1 ) « V TR AN I S) ) / 

O'F(S) * 0F(4) - R- 


-C(26) 

I ) ) •iOC 


ecu ) )/ 

CUO) )- 


HERE AMEN REAOY TO PRINT- 
SEE IF TTY mode. 
IFllTTY .Eo- .0) 60 TO AlO 

HERE IF TTY OUTPUT. 


Kt 1C TTY VUTHU I . 

WRITE (6.700) ex r »CCl ) ) ,CC(2 I ,cc( 3) »F (9 ) ,F ! 12 1 ,F ( 13) » 
b CH(4I,Ft7).CPa,C?T,FUl,FU7),?T02,VI,VE.CUl), 

6 FRE0.TVnT.AVO2DF,RHT(21 ,C( 10) 


700 F0RrtAT(7F9.4) 
return . 



40 ! 



fiu TO \??0 


610 IF (H •NC. Hi fiU TO \??0 

— N » 0 

naiFE l6tl80S) 

- >220 H « M + 1 
C 

-rt«|T£ (A,1810» CXTi H(5t »Ft5* 

C 

---JfRlTE (o.tatS) - IClll. I, ■ ltil» lOCtl), I o 1*3)* Ft71* F(l)* 

1 PA.12 

*<r>>7E (6,1820) C<.tl)f F(9). F(IO). KW), FtDi PAU2. tH.( 1) . 
t CPH( 1 ) , Cri3( 1 I 

iTRlTE (0,1323) (C(n. 1 * t, 6), (OCIilf ! ” H,6)* CPd* Fll7), 

I F ( 1 d 1 . CH( 2 ) , CPH (2 ) 

*!|R 1 TE (o , 1830) ICtl), 1 * 7,V). (t)C(ll* I “ 7,9), CPT, pTo2, 

I PTf-2 , CH ( 3 ) . CPh ( 3 ) 

_ flRjTE ( 6,1833) (OCd), I ■ 12<14), (C(I), I » 12,18), CH(*(), 

I PHCSF 

WRITE (6,1880) cc(2), F(12), C<6), CPu, F(17), F(lo), riVB^ 

1 PkvB, CHdC2) 

WR I 7E- ( 6 , 1 88S } CC(3), F(13), C ( 9 ) i Cl'T, PT02, PTH2, HVT, - r 

I PHVT, CHa(3) 

(,rt!TS (o.iaso) (T.AUd), J l.S), VI, VE, Cdol, dll), DCUO)*- 

1 oc t I 1 ) 

..RITE (6,1863) FRi:Cl.TvtJT.f>EAOVT,HRATC.AV020F,nSVOL 
1230 KETURM 

- - 1270 FOfit'SAT (SH XXXX6A7F1.0, 8 > 

1292 Format i < jf i o«8 ) 

• -- 1806 FORMAT H HI > 

1810 format 1 1 fiOoA'8HT 1 mEF 1 M , 8 , 78XoHAi.V R«F I C . 8 , 3 A7 HR tl 01FF,F3.87 
_ I 1 6A3hC02ax2n026x2H-j27X2 1 riO t R 1 V A T I V £ STX8HPC026X 

2 3Hp627A3riPtJ27A8HlH+)7X2HP),6X8HHE)02) 

ISIS FORMAT ( 3X5H-1, veai.AK9F 1 0, 8) 

1820 FORMAT 1 SABHARTtRl Ai,3if I O. 8 , 30X t6F I a, 8 , F6 . 8 ) 

1S2S FORMAT ( AXSwOftA : M 1 1 F I n. 8 ) . 

1830 format [ SaAMT 1 SSbE IIF I 0.8 ) 

IS3S Format (dX3HCSF3QA8F lf.»8l 

1S80 FORMAT (8a7mV wR A I HOf'] C . 8 . 30 X . 5 F 1 D • 8 , f « • 8 1 

1886 format (3X8hT T. I S6U£3 , 1 0 . 8 . 3,0 X t SF 1 0 , 8 , F t) • 8 ) 

laSO FORMAT 5 5X 1 BMTR/.I.SPORT TIP'S 8 X 7 h A bt> A 2 H V b b X 2H V T8 a 2 H A T a )(Z H AC2 X 
.. 1 2M« ♦8 ,x2mV I 8X2HVE8A 1 wQ9X?HFb7X 1 1 hOEK 1 VhT 1 VES/Zl A , I'OFl Q. 8 1 Fa , 8 ) 

135b Format ( 3X , 9HR3SP f R t.w .re. m,' 2 x , l 3 H!i.i.s:VTE Volume, fb,8, 

1 2A.8H0 S V£mT..F!i. 8 ,7X ,'l rirtP'tAKT KATE.f^d* 8, 

2 2X,7HAV020F,F8.8,2A,bhL(';vCiL,F&,8) 

-- C BATCH p.OOe ,-CRX CAkO REaO... 

c 

C AIlL ySF I.ORK CARO 1,1 th TlMi»C aS 1 .':l) 1 C nT 1 0 M 

c OF EnO of RuM otCA'USE 1 106 HAS PKObUtM 

C . . T. I TH EMn» or, READ, 

800 P.EAO { 6.,30p . E..0«’2 ) ViQ.K X 2 , i5UR A T 

300. FORMAT t F6 .2 .3 X ,F6.2 ) 

c 

IFIOURAT ,ST« C.) Go TO 6Q6 

C here if REaO IkOICATIO), of £M0 LF RL', II,' bATCH HOOt. 

.... Cl 15) a D, 

GO TO 1210 


606 •..RITE ( 6.3061 AO R< 2 , u„» A T ,C xT 




228 3C5 format ( *0’ * H3 (*••) / 

229 — 1 ‘ *CR< load CrtG. { « ,F A.2, • .vATTS FOK*. 

230 ■ 2 Fb.2,*HH'<SI at* .f 9.‘» . ’fUflS* ) 

231 i07 T I MEDFanURAT^CXT 

232 tImSDm«CXT 

• 233 C- SViTCri RFSPO><SF.Ss T 1 h£ Cm-.STA‘*TS Fo« >.*ORK LOADS A'^O TISSUE. 02 

23A C KETaoOLIC RaT£. 

235 • 1 F I «0riK2«£jt •t’ORK >RMTtJ2«Rh«T( 2 » 

235 C DECfiEASJ'.G lORK LOADS. - 

237' - • Ir(A0RK2.uT.urjRK) Ri-iTH = rlriT(2) 

238 IF{.i0«K2.LT,.-.0r'(;)KI\Tts = SSU2‘<tA0RK2) 

23V lFt(wORi<2.LT.AUI(rO.ANn.(.JOR‘<"G£*&0.>> TCT‘=2.3/t2.»y.ORK7200.) 

2*iO |F(tA 0RK2.LT.A aKA)»ANti.(.-(0R>;.LT.S0»I>TCT‘«‘).5 

2ai I F ( 7 • . WORK ) GOTO! - -- 

2A2 ■ tOl .iOR<a..OS»:2 

2h 3 - MARK£R«n 

2M‘4 (..•R£b7«2 

-• - 2M5 "C Tissue 02 metabolic rate. 

215 RMTt2)«RrtT8«(RrtTB-RMTM)*eXPl“TCT»(CXT~TlMEOfJ)**50) 

2H7 - VTIME" 1 . 1-1 . 1 .Ear ( -TCt« < CST-T IHEON) /3 .85 I • 

258 C TERM USED l;. VI TH»T IS A COHPOliEMT OF TKAI.SIENT RESPOtlSE RELATED 

■■ • 2-.V C TO WORK LOAO. - 

• 2b0 R’ILT. =RriT8- iRMTii-RHTKi • ( 1 .-VTIMU) 

■ -26.1- - [F’VTiME.&E.l.) R.ML1< = HUt8 

252 C TISSUE CC2 METABOLIC SATf. 

253' RMT ( 1 ) a . 83 »rmT I 2 ) • 

255 1FITV-IT.GT.3/. ) RKT ( I ) = t T V.MT * 50 . ? 7 > • R « T ( 2 ) / 8 8 . S 

255 IF {Cl3S).LT.C(50M G0T02 ~ ■ 

255 aRITE lo.333J RHTl U IRKTI2I 

257 — ^5X02 

258 End 





0B6-G03H32*TPPS.SS02rt 
j FonCTJON S50Za(aJ 

2 C CALCOUATION of S7EAl)T-3TATe OaY<jE''I REOU I RtiSENlb FOK VAKIOOS LEVELS 

3 C -- OF nOR(< LOAD (X“wA7TS). 

4 Si COHMOi.'/RlNTR/ROljT U 0) , C I N ( 1 0 J 

5 VOzROT-C IN < 3) 

6 S3o2.v = V02R0T-.0£00+<.n00A'\5031&»o.l<;«X(/.2S 

7 — rlETUKH — 

8 ENO 


CPRTiS terg 



OBi-GOiMSZ^TPrs.TERG 

J , . SUHaOUTl^E T£«G , , 

2 C»»* GE cardiovascular UCNP model 10/73/73 

3 coMHON/sTATe/x<50) .xonrtsoi 

R 2/STATE/LRA.oRu,r;LAHiLV.3PA.OPCtOPV(0AA.0AI?C»C!LAAiUUTAnilLTA»DUA0Af 

5 30lAOA.OClLU.<)LGSA.CLGAR,fjLGCAP,i;LiVE,DUGSVrQFey,oAiiVC,yTnVC,aSPVC, 

i ‘tULOC»UllPc*OHCAP.ynSv,OJV,f5cORiL'cbMA.C'I-‘-A.Lc'3llV.UPOVi>VlMV, 

7 5‘i''£^A .OfiALE , *AE^ V ,1«Et ,0D < 10) , C»SAS 

6 P/S'ATE/CRAtCRv.tLAfCLV.CPA.CPCiCPViCAAiCARCiCLAA.CUTAjCLrAiCUABA* 

— 9 7CLAi!A ,Cf 1 LL .CUUSA ,ClGA« .CL6VE .ClGSV.CFEV ,CABVC iCThVC tCapvC * - 

10 acLnc.cupc.c’fSv,cjj,cc«MV,ci'''V,cpOV( 

11 9CRE lA .C~£NV ,Ci>< lol 

12 A/SrATE/pRA.pRv,f-LAiPuu.PPA,PPCiPPV.PAA,PARC,PLAAtPUTA,PuTAipUAbA, 

-13 . faPLAbA,PclLL,PLGS)A,PL<3sR,PLGVE,Pu«SV,pPtv (PABVCtPTHVC.PbPVC* - - • 

I*) CPLOC >FUPC »PiSV ,PJV tPCF''! / , P i llV , PpOV . 

15 — OPH£N A »PRE‘i V ,Po l I 6 I , , p.i_c " • 

li COMKOD/qT ATE/ , 

17 - ERPA.RPV.RHV.riAVtRPA.KpC.PPV.RARCiRLAAiRUTfl.ftUTA.RUAeAt 

18 _ PRLAtlA.RClLL.RuGSA.f5uGAR.RLGCAp,RLGVE,RLGSV.HFEV,KABv'C, 

_ . ] , -GRThVC ,RSPVC ,Ht ()C,RUPC,RKCAP>h«SV,RJV,RCUR,fiCSHA,ftiHA»RCSMV, 

20 hRPoV,RlMV,PRCfJA,rtrALt.RREFFiRRt;NV,RD(ll).RS<& 

21 - i /ST A TE/FuPA ,Fla A .PLAkC .FLLA A , Ft,i;T A ,FLLT A I FUUABA , 

22 viFLl. AfaA,FLCILls,FLC5MA,rt.!HA,FLRElvA,FLDI'!(bl 

23 K/ST aTE/ vl 50 ) 1 VU 150 I ,Pv- I 3*1) . PEXT I 32 > I E ( A 1 -■ ■ — 

79 •.PRMiABlAS,TB[A5.TThA7.TM0')EL.Sp«CE(5) 

25 L ,2 (90) ,«< ( 201 ,hR,5V,C0 ,RT ,PtX ,t.,PSYS ,PUYS ,FR£e - — 

24 K.(,02i/OT,AVD,PlAO.PlTM,?rtP.THETA,SF 

27 - N.TT0TiTAS.TyS»Cl,C2)&yE.(,P£XlN,Tft 

2B • lOUMKY ( I 3) .T.DPHT . VLEG 

29 - CALL Ain - 

30 Call Cc;jTRl 

31— — -call CVS — • - ' - 

32 1 CALL ALGO(T) 

- 33 Jp (T.Gt.TTHAZ) TmETAbO. ■ 

39 IF (T,LT.AiC{ 20) 1 GO To 1 

35 —■ “ ■ call exec - 

36 RETURN 











08i-&03M32»TPFS .CVS 

j bUftROUTlNE CV5 

2 C 6E CArtOluVASrOUAR LBljP t'llJt.'t.L 

- 3 — C Cu>iTI!ClI-EU SYSTfH 

*i C0MH0»/RINT«/'Rn-tl01iPaijr(!0i 

5 - • COMHOt)/S’*TE/ X t bC ) , AO/^T ( ST ) 

6 ^/STATL/ORA,'J■^V,.JC.^,.l,v,lJP/. .•»PC,iiPV,.»AA,UA»<C.WI.AA,QyTA,'<»LTA,ijUAl>A, 

— '■ — 7— 3 CL Ab A . CC I UL > vL 6b A • <«LS aR . 'iiL t'C A P < SLC V t > U b V » OP t' V « 0 A tl V C > xT H V C I Q 5P VC • 

8 HOLOC.OiUPC«0><CAP*JnSv,pJv,'3COH,OcS.''A,HI'tA,.7cSf1V«'JPOV,i,iI«v, 

“ ^ ' StfREHAiQPAL£t»-»iEnVfORrT*»OD( !Q^ 

tn 6/STAT£/CRAtC<V.CLA.CuV.CPA.CPC.CPv,CAA,LARC.CLAA*CUTHiCLTA,CUAflAt 

U • 7CLAo-.CC 1 LL .Cu»S^ .CLCaK ,CLGVE ,Cl bSV .err V .CAAVcCTH VC .CSH VC I 

!2 3CL0C .Core . CnS V , e J V .CCSHV ,C t pV ,CPUv I 

13 9CRe‘JA.C«ENv,C0( irt» 

I*) A/STATe/PRA,PRV.PLA*P(.V.PPA.PPCiPPV,PAA,PAi<C>PLAA,P<JTA.PLTA*pUAoAi 

is — 8 PLAUA. PC I LL .PLLSa ,PL6aR .PL&V t ,PLbSV , PFLV , PAoVC . PTHVC (PbPVCi 

J i CPLOC.PlJPC .PHbV .P JV .PCS'tV ,P I PV ,f‘pO V , 

17 OPKEHA.P,<t-tV ,Pfi 1 I 6 1 «Hm ,P«c 

la co.AHo:*/sr ATE/ 

19“ ERRa . KHV ,r(HV , r! A V » «PA . KpC , RPV ,K arc * KL AA » rut A ( RLT A t riUAb A I — - 

20 FKLA3A.RClLL.RL6&A,riLGAR,K|.C.CAp,7L'>VfctUL6SV.RrEV,f<ABVC, 

21 • ■ •- GrtTiHVC ,«SP VC .P 1 .OC ,nUPC .HriCAP .RriSv , P J w , RLuH .RCSPA iR INA I RCSHV * - 

22 HKP0V,Rli'lV,R(«e.iA,'<RALL,R.<FFr,Kr<EuV,rtJ(tll>K5»;B 

23 I /STaT£/fLP A. r L aA iFLAKC .FLLAA, FL uTA , F LLT A .FLOaB A » - 

2*< JFLLA3A.FLClLL.FLCbHA.FLlMA.FLKEtiAiFLDHl3) 

2S K/StAtE/V t 50 1 .VL ! b-i I .Pp. ( 3't ) .PE/ 1 ( 32 1 ,E ( A) . _ . 

28 ..PRJi.ARIAS.TaiAS.TTriAT.TfnoEL.SPACECbJ 

- -27 - -- U . 2l **3> ,v,K < 20> , SV , C.T.RT. PEX , , PS TS tPuVS .F«EO - - 

23 M.V025QT.AV0,PlAb,?lTn,PMP,TbtTA,5F 

29 N*TT0T»TAS,TVS,Cl.C2.&ME«,PEXifl,TR 

30 • .OU.TMY ( 13 ) .T .OPKT . VLEr. 

31 — OlotRSIOM PR3 ( I > ,CHM i2 J ,P&0( 50 » iF IinR ( I 2 ) - • - 

32 E;u 1VALEI*CE lPSS,PSA).lc<Pll>.CRrt).lK50!l).«RA),(FlN«in .FLPA 1 

33- -- i, , tPU{ 3 I . T T > . ( PO t H ) , TSVF. I , 1 PU (5 1 .TRSP ) . I POl 6 > ,ThP) 

3A 8 .(P0<71jTPS) , IPOiai .p2) 

- - 3S . c T IS euaPSEO time - •• • 

36 C TT IS A CLOCv; FOS ONE bEaT 

37 - TT.T-TSVC 

38 IF ITT-TTOT) 1002,1001.1001 

■ 39 - lOOl- TSV£=«T 

hq C*«** 

- •- 'll - -- 0»M0« 1 7oQ • + ‘4bnO • /3 . . 1 VU200T-. 5 ) ■ , ■ ■ 

H2 0?£R = 3PM.>-P0 110) 

- - N3'— 1 F J OPE=? .LT .-bo . 1 PO 1 1 1 1 bPo t 1 1 t --OO 1 • - ■ - 

PR IF I0PER.3T .50, )P0 1 11 ) =Pr) t I 1 )*.00l 

— R5 • " - I F [ V0200T.LT . .b I POt I 1 I »0. 

M6 CO«X { 33 1 /TTOT • ,06 

- - R7 /( 33) oO.O 

RR Pn«A I 1 0 ;/TTO f 

— R 9 — - - — — - • X(|0)^t3,0 - — — — - .. 

30 PflC'Xl 1 3> /TTOT 

51 - • X( 13)-0*0 

; 52 P0( I )“*13R)/TT0T 

53 Xt3R)“0.0 

SR SV»TTaT/60.»cb 

56. RT»?3t 1 )/CC 

56 0 5 FF“-V t 50 ) *•/ ( RVl 




C" 26 
26 


XI iB)«x( iej+oiFr»o.6 
X< tV>»X< 191+DIFF'O*** • 

PSyS»SyS 

■ PorscoY<? 

CA1.L X!0 
POt 10) a C ♦ 

CALU COmTRU 

-TE'iPcTenP+0‘2 — 

IF (TEMP-TJ UO.lU.lTl 

COfjT 1 NUE 

CAt,L EXEC 

SY<3»Ci«0 

OYS^IOCO* 

TTq 7«60./HR 
T A5"0 . 1 0+C.O?«TTOT 
TVSaJ.) 6+0.20»T70T 

IF IT.LTi**!. .OR* T.GT.'IJ.) GO TO 20 
Oo 10 J=1.32 ' 

PGn)*=SIhtTHETA/S7.29«lB)*Z(I>»1.05«980./1332. 

COnTUjUE 

IF { TMnnSL.GT.O. 1 • 60 TO 26 - 

GO TO 30 

IF 1 A3S 1 THETA » .UT . 1 .E-5 .0«. T.lT.'JQ.I GO T030 
IF UbSlTHETA) .ST. l.E-SI TILTD»1. 

IF I A3<> ( theta ) ,6T , I ,E-6) GO T0"30 

IF (T ILT0.GT.2. 1 GO TO 30 

00 Zo jsl .32 

PGU )=0. 

TlUTOas. 

CONT I tIUE 

VUEG=0* 

DO 201 1=15,20 
VLElj^yLE^i^Vl I ) 

COnTIMjF. 

VLEG^vlEG^VUI 1 e !-VUt 1 9)-VU I 20) 

T£mPV=0. 

00 1 6 1 = I » 32 - ' - • ■ 

T£KPV=TeMPV+/U( 1 > 

SPACEl3)=V(50)-VL,E6r*Tp«PVTVU!ie)*VUn9)TVDI20) 

i-*VU(13)-*VU(16)'*VOlI7) 

COnTK.OF 
!FtTr-TAS) I ,2,2 

SAS = S JK( 3. 1 9 I ft'^TT/TAS >• • - - 

£ t 1 ) =0.05*0.wS.SAS»SF 
E(3I°C. 1 2.-0. m»bAS*5F 
KSpVC‘=l20.*SAS»*lO. 1/1332. 
r<TrtVC«t 10,»SAS»20. 1/133 2. 

GQ TO 3 
e 1 1 1 “0. 05 
E ( 3) =0. l 2 
RSPVC«.0 1 501b 
KTrtVC«. 0075075 
TV»TT-0. 1 

IF ITv'.LT.O.O)TV = 0.0 

iF<Tv«T\/S(5,5,5 

SVS"SIM!3.1‘U6'*TV/TVS) 





ss- 


e(ZJ=*Q.ni/S + .3V*Sf«!,vs 

E(H}*0,O2*l.i>0«SF»SVS 

^0 TO & 

e{2>"0.OJ75 

E(‘n»6,r)2 

COMT I -'(•JE 

00 11 I = l .<< 

CrtPl I >•!./£ H > 

iFixiH) .Lr.o.onu*u»o,o , 

COMPOTE VOLUMES 

V (SO)aO.O 
00 5S l«li32 
V( JJaVJll l+XC I ) 

V ( so ) » •; < so 1 *v { I I 

V(SOMvl5Ql-V(9>-vllll-Vll3>“VtlO)-vUn«5l-VUU9)-VU(20l 
ir ( fn£TA»ST»Mt>iAL'0*T*<jT»M0*) PT TM^"2*5 
MESPIKATOKT 'oUMPS 
IF {PEX.E3.0.0J uO TO US 

) SO TO I IS 


I3Z 



lFlT.i,E.40.a •Oil. TmETA.lT*' 

133 - 


— 

-T«5P = T!?5?*T-TP5 

13M 



lF{T.?SP.ST.TH»TaSP»3.fl 

135 


- 

T1»T45?/T« 

136 



PlTK“-2*6/*“l9.’7a4»Tl+56»409' 

137 



- 

OE?Tri»< V02i)0T-l .1/2. 

1 3S 



IF(DEOTH.LT,0.l0EPTri“f), 

- 13?" 

— 

— 

lFU)£PTM.5T.l.5)OEPTri*loS 

IlO 



P I Thap I TH-OEPTH 

1 M ! 



PlA8“-PlTH/2« 

1M2 


1 IS 

C04T 1 f.OE 

.143 

* 


00 7 1 iB 1 . 1 2 

144 


7l 

PExTt I laPITH 

14S- 



— 

PEXT(221“P1TH 

146 



PEXT (23 1 “PlTrt 

- -147 

■ - 


00 72 Ia2rf»32 

148 


72 

PEXT 1 I I aP 1 AS 

1 4? 

- 


PEXT ( 1 4 ) a? 1 A6 • 

ISO 



PeXT ( 2 1.) aP I Ad 

■- 151 

' *C- 

— 

MUSClE pu«p 

IS2 



THP^THP-fT-TPS ' , 

•' 153 

— 

“ 

tPS“T 

154 



IF (TMP.ae.l.l TilP*0* 

155 

- 

- 

SttPaS 1 (l.( Z * »3» 1 4 r6» TrtP j 

is6 



P M ? ■ 4 0 . • S .4 P 

— 157 • 

— 

— 

IF<TrfETA.l,T. IS. )p H pa 1P.*SMP 

laa 



IF (SmP.lT.0.1 PHP“0. 

- 159 



IF (PcX.lT.1«I PM?“0« 

160 



DO 44 lal 6, 19 

lot 


• HH 

PEXTl I lapMP 
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COMPUTE PRESSURES 

■--—■16 3 -• 


. - 

Pl»P2 

164 



P2-PL V 

165 



00 12 I-1.7 

1 6 


12 

PRSli)«xll)/CMP(I)*PExTin 

167 


... . 

0P0T“1PLV-PI>/12. • .0Q> ) 

ieS 



IF (CPOT.ST.PO ( te n PO ( 1 0) aOP! 

- 169 

— 

» - 

00 13 I a 1 5 . 1 7 

170 
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PT IS*SPACE i ‘I J 
— PSBlASaSPACElB) • 

00 IS t*l8,20 

PRSU>»Ani/’''0(l)«2, + PEXTn>+PTlS*PGBIAi 
5 IFU ( n .gT. VU ( I ) IPRS I I ) B 
& (X(l)-yU<ll)/CHPIlJ-*.PEXT(I) *PT IS + PGtJU 
DO iR j.zR.oa 

^ PRSi I > “X 1 I > Xc^PU <*PeKT< 1 > 

PAA*Xt9»/CAA + PI7ti 
PUTA»Xl?9>XCUTA*HITH 
PLTABXt ^^J/CLTA■^PtTH 
lF(piJTA.GT.S’fSl!>YS"pUTA 
IFfPUTA.LT^OYSluYSBPUTA 

PLABA»PjAB-il.624*0»OP2?6S*V<lR)*0»0097( 
PtAbA«xj I‘t)/C(.A&A*PIAR 

• AtjCOP.lNAL VEf<A CAVA 

PA(WC““G.R99A + CI.Cd2AOA»yt21)-0*00033S98« 
.*O.OOCnoO‘»SC'26*V'(2J J«Vl21»»V<21 J 
IF ( A ( 2 1 1 • GT t 2C0 • (AND. X I 2 1 ) (Ut • 3S0« ) 

, PA5VC«:.3*</l50.»lXt21J-?00«> + ltl5 

• Thoracic vesja cava 

PTnyC = -S.b006»0*l 15‘t*v(22)“0*000<>587a*V( 
,+0.O0Coni236»Vl22)*Vl?2>»VU2) 

• IF (a( 22J*GT. IbO. .ANpi. X ( 2 2 ) .i T . 25 0 • ) 

» ?THVC = . 3/ 1 00 . « tX(2?:-15D«> + UJ6 

- PSPVC>-3.R99 9>0.92'TOV.X(23)-ai0HZ2‘l6«X<: 
,to.OSO«j3i&5«a (2iJ*Xl23J«X(23» 

PT« VC = p TH VC *PEXT t 22> TTS I A5 ■■ ' 

PABVC=PAaVC+PEXT(21)*AaiAS 

PSPVC = P'5PVC-»PEATt23> 

QRA"(PPA-PRV) /raa 

FiEAKT t'.OOEL - • ■ 

IFCPRA.I T»PR* IDRAaD.O 
QRV«=X ( 09 I /FlPA 
1f1SRV,lT*0.0)«KV=0.0 
XOOT t09 ) =PR .'“?PA»KRV*qRV ■ 

JFUOt)T<09).l.T.O.O.AIin*ORV.EQ.O.O) XDOT tC 
•' GUA“<PUA*Pl-Vl /Ri;v 

1F<PLA.(.T.PlV)0LA»0,0 
QUVBX ( I 1 ) /FUAA 
lFlOUV.LT«n.O)fSl.VBO.O 
-- XDOTtllJsPl V-PAA + PCtdi-KAV^ClLV 

Ip{X£jOT( I 1 > •i.T.O.'’.AfJn«OLV.EU.O.OIXoOTt 1 
-- • -fuLonsAHY Circulation 

^ 5PA3 tPP A-PPC ) /RPA 
0PC« iPPc-PPv) /KPC 
0?y» (PPV-PLA) /HPV 

arterial iiOOru 

GAA3<PaA-PuT^*Pv112))/HuTA 
■■ U JT A? < PUT A-PLTA»PG < 1 3 ) J /RLT A 
eLTA=tP(.TA-Pi.A3A + pu( In I J/RLABA 
0LAaAa(PLA5A”PC!LL*Putlft))/RClLL 
LEGS 

QClLL«(PCILL-^PGtl6)-Pir.SAl/RL65A 
Oi,GSA» ( puSSA-plGAR ) /rtt Gar 
--■•OlGCAPb tPLGAR-PcGvE) / rlGCAP 
RLGVEa.OS 




228 

229 

230 

.231 - 

232 

233 
23** 

-23S — 
23& 
'237 

238 

239 • 
2**0 

241 

242 
2*f3 
244 
Z4S 
246 

' 247 ■ - 

243 

249 

250 

251 - 

252 

253 

254 

255 

256 

257 ~ 

258 

259 — 

260 
261 - 
262 

263 

264 

265 - - 


IF tgLG\/E*LT«J«0)KL<9«/L>67. 56/56/ 

0LGV£«( PLSVE-PLJiStf ) /Kl 6VC 
RLSSVa.OS 

IP [0LGSV.UT.0«0iHUG5V«67. 567567 
iSLSSV^t PlSS V-P& (l9)-f'FEV)/Rl-GSv' 

9e4-JuS MODEL 

8FEV«.r721 

IFOFE-/.LT*-''*C)RFc.v»67. 56 7567 

QFEv=(?r£v-PG(2oi-PA&vC)/Rrev 

— aA9VC”lPA2VC-PGl2l l-PTHvc»/8AaVC 
QTHVc<*(PTH/C«?G(22)-PRA)/KThVC 
GSPVC«>.<PSPVC-P<j<23I-PRA)/RSPVC 

HE A Jt- ARMS 

OLOCa (PA A*PS* 24 J-PLOC ) /RLOC 

iJ?C»(PLQC-POPCi/RUPC 

u6RAIM»SIS(3t* 1000. /6n. 

UAR40 17 .25 

— OHCAP-'SdRA I «*(2ARM 

QhC AP» ( PUPC-PH3/ ) /ftrlCAP 

ariSV»(pMSV-pJV)/RnS/ 

RJVx. 004301 

1 F I <JJV.L7.0.3> HJi/»67.567567 

0J7*(PJV»?G(27)“PSpvCl/RJV 

CORONAH/ CIRCDLATlOt. 
OCOR=.(PAA-PPaJ/KCCR 

COMriNUlTY Fnfl YFMOOS RETURN 

ORET»GSPVC*jTh»C*.CDK 

HtPAT ic-sPLAMCrtNic circulation 

0CSMA!>(PL7A~?CS!;V ) /SL5 *A 
- U(CSMV«<PCS4V-PP0/ ) /RC5HV 
QP0V»CPPOV-prHV’C>/RPOV 

KENAL CIhCi^i-ATIOM 
QRENA«lPLA6A-PR£f<A>/HRENA 
' QRALE*< PREnA-PREnV I / (R fiALE*RREFF I 
aRENV={ PAENV-P AfcVC ) / hREUV 

— SAEuTOM.pOnE 4ARR0«VtANt> OTHER 
0S<aa(PLAoA-PA6VC)/RSKa 

■ • state VAKinoLt derivatives 


266 

XOOT ( 1 1 ajHE7-Ci«A 

- — 267 

XOOT ( 2 1 •3RA--RV 

26S 

XOOT ( 3 I aSPV-OLA 

264 

XOOT ( 4 1 k9LA-'3l V 

270 

XOOT <5 Jn'jrtV-.PA 

27 1 

X00T(6) = ‘.P4~mPC 

272 

XOOT ( 7 ) SwPC-vPV 

273 

■ -- XOOT (t J=uUV-wAA- 

274 

XOOT ( 10 1 aPAA 


275 - 

276 

2 7 7 - — 
273 
279 
260 
23 1 
232 

-233 

284 


XOOT( lZ>BLoTA-aLTA-0C5HA 

6D0T(14 1».LTA-',L«»A-i,»ENA-f;SK8 

XDOTl i5}6v,LAoA“.;C(LL 

AQuT ! 1 6 1 = JC: wL->.u6Sa 

XOQT ( 1 7 lawLGSA-.'XLoCAP 

AOOT I 1 ft > 6 NL.CiCAP“CL6V£ 

XOOT ( i 9 l=aLfwF -wLt.SV 

XDOT ( 20 1 aiLGSV- ;r t V 

XOOT (21 >PwFEV*w*.i>v'C.'.AtNV*«.SAA 

XDOT(22)6»Ag/C*OPCV-CTHVC 




zas 

XOOT { 231-ajV-OSPVC 

286- 

XOOT (29 loiLOC-fOPC • 

237 

XD0T{ Z5J »3UPC-0mCAP 

2 3 8 ■ 

XCOT ( 26 ) »4hcaP'’OHSV 

2S9 

XOOT(Z7)=«ihSV“'JJV 

270 - •- 

XOCT < 23 ) ■’.JcS«A"f>CSMV 

271 

XOqT 129) “iAA-Quia 

292 - -- 

— XOOTt30) = JC5*tV-OPOV- 

273 

XD0T{3l>»3ReNA-uRAl.e 

27H — • 

XOOT (32 )»'.'9ALE-wfiENV 

295 

XDOT (33 ) «3L V 

296 • — 

X00T( 131-PLOC 

297 

XOOT(39)»?UTA 


^ iSCTIlOU . -* . 

* 9 g - ■ 

299 

nt 1 wrtJTf . 

£ND 


IPST.S CONTSL 






o 



0a6*G03'l32»TPr»«C0MRl. 

j- SOBROUT|(ve CONTkL 

2 C0KH0N/STAT£/X t 50J » aOoT ( SI), 

3 . . 2 /STfcT£/'uSA*'i'<V,w!,Ac}t\/iiJPAiUPCi'JpV,.jAA»tfAMC if»LA**wOTA,wi.TA*i(UAaA, 

*1 3uLASA,TClLL.Jt.S!»A,3LGAR.^l.fiCAf',0tM/E.*iJI-uSV,UFeV,uAc3(/CiuiTHVC,0SPVC* 

- • 5 — HOUOC * JiJPC .OBCAP .iJnSv ,OJ V ,5 CoH ,HC>IIA I "A »0C5HV f t4KOV I -( I'ltV , 

6 50i?fC>4At'3RAI.Ii*'3'<E!»''.0‘iLTiJl}tlCI,ISi<r' 

—7 • - 6/STATE/CAa ,CKV,Cun,CCW.CPA ,CHC *CP^ »CAA ,CAHC »CLAA ,CUTA ,CU-TA*Ct»AoA, 

8 7 CLA 3 AiCCtH.,CLGSA,CLG 4 ^,CI.<'w£,CL‘>SV,CFEV,CAiivC»CTnVCfCSPVC* 

---• -- 9 - 8CUOC ,cjpc .ciSy ,c jv .ccqfi / . c I M'/ .cpov , 

10 9CBE.’iA.Crt£ .V.COi Idl 

-II A/STaTE/PB A « V » PLV 1 PPA , PPC I PP '' t P A A ,PASC ►PLAh , PJ r A ,PUT A * PUaB A » 

12 b?l.ABA,PCIv.L,Pl5SA,PLSA«.PL'iVh,PLUSV,PFCViPAavCtPrHVCtPS>PVCi 

13 CPLoc .popc tPHS<'iPjy I PI MV ,ppov , 

m 0P«E'4A,P«ESV,P0l lA) .Pi.PmC 

15 -- - , COHilON/ST AT£/ 

16 . £.IBa .SRV ♦ =» IV t BaV . «PA , ><PC • 'IPV ,R A'lC > ‘IEa A * MtlT A ii'<LT A » hOABrt • ' 

J7 FBLA3A ,RC ILl t«^U63A ifiLGAB .RLCiCAP irtL'iVEi Rl.GSV tRPE V iKAoVC , 

18 GBTriVC,R8PVC.'<LqC>-(JPC,BHCAP»KH5y*f<JV(KCl>K»r'fCSMA»f<l«AiRCSliVt 

)9 HRPOV iBtMV ,RPu»jA ,BBAut ,'iH£FF .RrKhV (Ri>( i I ) »RSSB 

20 l/STATE/FLPA.Fl,AA.FUABr.FLLAA.Fl.uTAirLUTAiFUijA0At 

21 - - JPULASA ,Fl.C ILU .FwCS-lA . FH MA , FLREfiA .rua.'l I d 1 

22 ' iC/STATfc/ V ISO) . VLt'lGOI ,?r. I 3Ml . RtXT ( 32 ) .t i 4 J 

_ _ , 23 • iPR.J , Art 1 A5 • Td IAS 1 T T«A7 . THOOEU .SPACE (SI 

2rt , L,Z ItO J ..,1(120 1 .Hri.SK .CoIkT.PEX .'.r.'PSYS, P07S»F«E0 

2S M. V0200T , A VO.P I Ad .? 1 Tn .PrtP . theta ,SF - 

26 «.TT07,TaS,TVS.CI,C2.6mEA .PEXINi'TK 

27 - •.DUMIIYI I 3> .T.OPBT.VLEr 

28 CC CVS-HE5P. INTEKFACE 

29 C - dLOCK OATA FOR IhTERFaCE |A 6 OuT 

30 COHrtOi(/R ! NTR/RINU 0 I I ROim 10 > 

. 3i_ . _ c INPUT FRi).N RESP. 

32 V02dOTaRIfi| I J 

• - -33” — "-’FRE«,"*R 1 n( 2 > ^ 

34 C 28 RA I N3«lf,( 3 I •1000./6n> 

--.33 - - - PCOZorUNl'U ' ■ -- 

36 POZ«RIII(SI 

37 ■ — 0PC0'2opC02-37.57Sl 

3s DP02'opo2nl07*9‘l82 

: 39 RVaCT» 1 .AUO?C02-rOP02l«0. 01273) 

.40 C ' RVAC\«l. 

, Ai - C- OUTPUT TO RESP» 

• 42 ROUT ( i > “CO 

A3 rout t 2 1 »>V0200T . _ 

44 rtOuTl3loREST02 

- - -- H5 - C ■ CVS HOOEU 

46 REAL Nurt9l)I tNUM9d 

47 . - 63UlVAL£:.CE 

48 2 ( ACCHET ,X ( 4 1 )) . UNh . X ( 43) 1 » I DA . X I 44 ) ) » ( DL »x I 4Sn . 

49 — 3 IDO 2 .X (4Sn.lX!.3.Xl49|).<POl2l,FLA8},tPl,l9).OTSl.lPO(l2I.OMS) 

50 COMMON/OtLA YC/ AVdTS ISO >. V 02 Ts ! 50 > t SAVE t 1 0 ) »F 1 t I S ) »F2 U.I 5 > *F3 I 45 I 

51 & .ANF ,T0N,F*tS.«£ST02 

52 C SAVE OtO XOOT (41-491 

53’ do 1 C 1-40, 4V 

54 10 5AVEU-39>oxU0TI 11 

55 -C' — 

'56 00 25 Ia| ,17.2 



HR^ iJAfE A6 -k/M|N 



S7 IF(T-»vK(m 26WSi2S 

SS 25 Wm*kU*1) . . . • 

59 24 IFU.> 27t27,28 

-• 60 • 27 PEXnO.O • - • 

41 GO TC 29 

■ 62 - -- 20 PEX»1 *0 

63 29 CoNTlriLKi 

. j -OXyGEN requirement function- V02U0T 

65 1 F I 002 . CT • 0 • 0 > 002>0 t 0 

46 -- VO2r<0T“.O00*<6BOel!J*ft7.2S 

67 PSft<=-l ,6 + 0Q2 

68 • £>T1«002 - .. 

69 OT 2” (2. •002-1. 27 51/1. 15 

70 DT3«0« + DL - - ... - 

71 OTtN«'S-.VlN(PS-v,ini .0T2) 

- — 72 - OT»FCN'SA(P£XiDT3iOT3.nTlN) 

73 C AUACTIC oxygen DEBT OA 

7** - - OAlh = . 15»(07-1 .51 + 1 .b .... 

75 DA 1 .= 5.1 1 N (PS V ,002 f OA IH) 

74 ; OAO«FCNS'MPEX ,0.0,0. ojOA 1 > •— 

77 T8«FCilS.VlPEX>0.Q.30Q.,2. ) 

73 — XDOT HRI-IDAO-UaJ/TB • — 

79 C LACTIC OXYGEN DEBT OL 

.. 30 OL 1 H = .8S* COT-1 «b > 

St OL I =3il I .N ( PS* ,0.0 ,OLlri) 

32 OLO«rCN5'«fPtA.O.O,O.OjOLl > — • - . — 

S3 T7 = r CigS.f t PEX ,C.C, 300« , 10. ) 

... -g,, xooT t<,5 1 => 1 olO-ou) /T a ■ ■ ■ 

65 c ARTERI aL-VENOUS OXYGEN DIFFERENCE AVD 

. .. 33 NUMVol«. 038*002 

37 CALL OELAYt0.O.O5,NuH9DI .AV0TS»'^UM90il ) 

33 A VD=*Y0200T/CU — - • 

S9 aDOT(4A>=*{nu;i 90-XM6) l/S. 

90 I r t PEX .£0 .0.0 1 FLA‘“*=0.n ■ • - . . . . . 

91 IF {FLAG.E J. I .0 J GO TO ftO 

92-_ 1F(PEX)60,60,41 • ■ ■■• ' • 

93 61 AnP*1.0 

95 FLA<5=»1.0 

■■ 96 60 I F 1 T .GT.TOU ) AmF-0.0 ... - . 

97 TAn»FC.NSm1ANF.3.,36*.3*J 

98 • XOOT IH9 )=• n 1 .CC*ANF-Xfj3 > /tan . .... 

99 XOOT l‘i3 > = t 5.5*PEX-XhH) /-6. 

lOI OM=002.25. /22. 

102 IF<PEX.GT.O.O)DTS = OT - - 

103 1 F I PCX .f, T . O . O 1 OnS=*Oi>l 

104 ... ■ IF 1 PtX .1 T . I . ) 0.'.«0.4S/0TS*DT 

105 CriC!lO.*J:in02/.0.3 

.... 106—.—— — 1 F ( C ME NO;. . G T >0 . S 1 C HEKdNso , 5 . . . ..... „ _ 

107 F.'* = 2 .•X.S'i*X‘iJ 

^ — JOB -- I F t FN.GT . 1 1 . ) rN“l 1 . •• ' - • 

’ 109 lF(P£X,GT..O. IF\S“FN 

no- IF (PEX.LT. 1 , )FN=Fi,S/0TS*OT 

til c controlled resistances 

1 1 2 ■■ -- c LEGS - 

113 RH£7ol50.-ACCH£T«S0. 




IIM IFt»?MET.LT. 

1 15 ROH»H50*-‘t5r'.»0(t/i>KllA 

116 IF (RDH.LT* IS. )nu!l3|5> 

1J7 (»L6CaP-(Kh£.T*R0m>/133?. 

113 RLfiARHoSIlET^KOM 

19- RLSAt(M«-Ft.»5H00./ U • 

20 RL<5AK» (rtLSAKh*KLGANH*=('(00. ) / I 33?. *R VACT 

2j C other eSANChES 

22 RHCAPO3670./1332. 

23 • - RCOH" I 205CC.-93V5.»0H/DMmX 1 / 1 332. 

2*1 RSK6>> ( 7640.-*-o000.«0L I / I 332* 

25. rtRALt* (5COO.- 1600. »DL) / 1 332 . 

2 6 O0MR»t)‘(/0KMX 

2 7 . - — jp (R'jhr.st . 1 . )Rortft*«i . 

28 RC5;iA«(5SOO.*IO/Q.«(Fij/n./2*+«oHR/2»J)/ 1332.»RVACT 

29 ' - - Son"0. 

30 00 92 t» 1 , 1 4 

31 FI ( 1 I“F1 1 1*U . 

32 92 SUH"F 1 ( t ) ♦SUH 

3 3 ^ PI I 15>>«CD 

34 COT= ( SUK + CO ) / IS. 

35 IF(C07.f.T.2S. )CCT"25. 

36 R?A = .0175-.007'S/21. «. COT 

37 RPC«.C59S“. 02*15/21. • COT 

36 KPV'RPA 

39 RREFF“0. - ' ‘ 

*10 IF (Theta. GT .15. .ai.u. t.rt.io.i go To ago 

m — . . (tCSHA= 1 2*1 ?0. ♦ I 477o.-*l I 30. ) • ( F;./ 1 l ♦ /2« •‘•tih/OHHX/2. n / I 332 . »RV aCT 

*.2 RLKAno (6000.-600. «PEa*RlGA9 N + RLG ARM 1/ 1332,*RVACT 

*<3 - -- KR ALE“ < 366o.'MflU2C.-S I 40. 1 •DL/C.9 ) / I 332. 

44 600 COnTII.OE 

45 SFa.6v*,374*(x(H0>-.9 ) 

46 tFISF.GT.l.l35)SF*=l.l3S».A6»lXl40)-2.143> 

47 — - • XOoT ( 40 > = 1.V02D0T-X { 40( I /T7 

4S 1 F (SF ,LT .0,67 )5Fo,67 

49 IF (SF.GT .3. >SF°3. 

50 IF (PEX.lt. 1. .AflC* ThETA.LT.MS. ) Sf«.48 

51 J F IThCTA .6T. IS. »'Ai;0. T ,lT .40. 1 SFn. M3 

52 SF«SF-.PD (11) 

53 610 CONTlmiE- . - 

54 C PRESSURE reference FUNCTION PR 

65 — PRN«66.*Cl •002+C2XACCNET 

56 EN»PRf>-PH/2.-Pf'C/2.-*XM3*J.«XN4 + FN/2. 

57 SuH=0 . 

56 00911-I.IM 

59 - -- — F2( I J BF2 U* 1 1 

60 91 SUM»F2 ( 1 ) +SUH 

61 Fansi^EN 

62 EP»(Si;m*En)/1S. 

63 P A VGsPO (11 

64 IF (PA\1G.LT»95* .AND. PAVG.GT.69.) GAlt.* 

65-6 ( 4 . I -3. 7 7/6 • • ( P A VG-oV . > I •GKt* 

66 IF (PAVG.tE . 69. ) G A I N = 4 . ! • GN £ W 

67 » - IF (PaVG.GT.96.) GA1I«»«33*C. 1. £A 

63 1F(PEA.G'T.O.)6AIgs.33.GNE.. 

69 ‘ - .D0PcC.533*(£k»Ga1/.) 

70 IF(DOP.lT.C.O)OUP**C.C 




J7l 



TOt«Ci000900P 

- • - 17.2' 

— • 

“ 

Hn»60./T0T 

173 

c 


COl.TROLLEO COMPLIANCES 

- 17 h 



tRC“lPK«-P(1Cl»*7 • 

*7S 



SOKoO . 

176 



00 90 I*!****! 

177 



F3 ( 1 ) «F3 1 I ♦ 1 > 

178 

— 

-90 

SUf-»f 3 t I > »SUM - - . - — 

1 79 



F3 I 1(5 ) BEKt 

160 

- 


e«C“ 1 SUM + ERC ) 

lai 



IF (ERC.LT.C.O J SO TO 7 

162 



IF<£t<C«GT.}i0.)EPC”60. 

133 



CL6v£«3.956*( 1«C-.0063*Erc) 

— 18*( 

— 


CLGSvb j . 1 HOS*! 1 •o-*oor3«ero 

16S 


7 

COnT 1 r,LE 

- - 16» 

- ■ • c 


respiration 

137 



ReST&2a.3l3 

-- J66 

- — 

- 

IF 1 THETA.GTt 13. .ANO, T.GT.80.1 REST02«>.37 

U9 

c 


FR£C"Vo20oT«6.2S*5.2a 

,93 

— 

- — 

IF CFKEC .GT. 30. 1 FRE0»3n, - - 

191 



TRb60./FRE0 

192 

• — 


AOoT (<* I 1 B ( vo2O0T-“0.38 }• .«(/3 00. 

193 



1F(PEX.EO.C.CUC>OT(‘<1|b>1./300. 

- ■' 1 9 It 

• 


;F(ACCMFT.LE.O.O.'ANO.P£X.ec.D.O)XOOT(9i)=0.0 

19S 

c 


oxvsen deficit function 

1 96 

. 


call DELAY<C.O,6,vO2OnT,V02TS,V02DD, n 

197 



X00Tt‘ie»Bi-YC.Z0L + PEX‘v02*DT + 0.33)/A0. 

196 



lF(002.LE.O.C.Ar.O.P£X.£&.0.0)XOoTM.'. >»0.0 

199 



00 31 InHC.S? 

200 


31 

Xtl )BXtI )*C.t*lXOOT( l}»SAVEn-39> 1 

201 



RETURN 

— 202 - 

• -• - 

— 

•END - - - • ■ 
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086-S03‘*3 2»TPFS. AL6C 


J 

2 

-- - -3 
H 

- ■ 5 
6 

7 

3 

• - 9 

10 
t i 
1 2 

I 3 

19 


- - iOpROUTlNE AU60<T> 

C iNTtOKATlON algorithm 

COmHON /STATE/ XtSO ,XDOT(SO) 

•4 l/lHENSIor< XOSCSC} 

Do 3 1-J .3** 

3 XO? ( H-XOOT C 11 

ti"0»C01 

1 Ft T.GT. ll . )Me.oC2 

T®T*n 

CALL CVS 
00 M r» 1 .3«i 

R X<i)"K/2.«<XD0T<I)*XDS(|ll*xU» 

RETURN — 

END 


OPRTiS XIO 



0B6-G03*t32*TPrs.X 10 

} SiUBKOUTlNE X 10 -- 

2 COMMO^/iTATE/X t6CD) 

• 3-- - cCKHON/'x I0D/N< 9) tNVfi 8) I IN 1 T • A< 9 ,6 > 

H *• DATA KT,.\i*Ti./'lHN.AH TIlT/ 

5 . _ T«X(59S J . 

6 IF tlNjT*Gl»0> GO TO ?00 

6 CAUL COATEIMO) 

• 9 call CTIMEthlTl ■ ■ . p . 

' 10 |.RITE(6 ,5JFD»KT 

11 G FORMAT!/' CARDIOVASCULAR TlLT EftGoMETRY MOD EL * 1 6 X » A 6 » * AT *,A6/- - 

12 • • REFER 70 GE-aGS USER CUloE TIR 79 I -MEO-* • • • » / / ) 

13 C - * '/TO simulate TILT EXPERIMENT ENTER 

19 XI 99b ) « 1 . 

— IS - C READ (S,6l X(99S) ■ — ■ 

lA C 6 F0R.9AT IrS.OI 

17 9 OR 1 TE t o . IC) - 

13 lO FORMAT ( >000 YOU 'A 1 SH TO CHANGE INITIALIZED DATA? <Y/n}») 

20 2q format I 1a1 ) 

— ■21 - IF IK.EO.KY) 60 TO 60 - • - -- 

,22 ARl’TE ( 6,30) 

23 -- 30 format ( 'OFLEASE ENTER !N0exU-600). VaLUEi CRt- <l3i£i2»6)‘) -- 

29 GO TO 90 

2S 3S r.RlTE 16,86) • -- - • • - - ■ ^ 

26 90 REAi; ( 5 ,&C . ESR"35 ) IiVAL.NE* 

27 -50 format I 13 ,£ 1 2.6* ' 

23 AHlTE l6,S5) IiVALnE*^ 

29 ■ S5 FORMAT 1 9X , 3H*« • . 1 9 ,F in. 9) - . - . 

30 IF lI.LT.l .OR, I .6T,600)'. GO TO 60 

3 1 XI 1) » VAUNE.Y - . 

• 32 60 TO 90 

• 33 60 ARITE ( 6,70) - 

39 7q foRpMat ( *0'jo You aish to modify the output li'stt iy/n>'3 

3S • read I 5,20) K . 

36 IF U.ec.KY) GO TO 20n 

•= 37 nRITS { 6,30) • — 

38 80 FORMAT COPLcASt ENTER POS 1 T I ON I 2-’ ) i I NOEX ( 1 t600 » , ' 

-39 ■ 6 • LA&El, CR; UI,I9.A6)*) 

90 GO TO 90 

91 as YpMITE <6,86) . -• 

92 86 format (' •READ ERROR*') . 

93 - 90 READ I5,iDG,EpR« 65) lP,i,N.VB ■ • — • — - - 

99 I QO format I I I , I 9 • A6 ) 

95. ARITE (6,1CD 1P.1.N.,R 

96 IQI FORMAT (9X,3h*** I 1 2« I 9, IX , A6) 

97 - - — IF tl'P.EB.O) GO TO 20n 

93 IF 1 Ip.EDt ) ) GO To 9 

99 — IF (IP.LT.2 ,0R. IP.GT.9) GO TO 55 . - _ 

SO IF I1,LT,1 ,UR* !,GT,A30) GO TO 85 

• 51 GO TO I9-3, lOZ, 103 . 109, 105. 106, lo7 , 108, 109) , IP 

52 102 Nt 1 l»I 

.. . S3 ... . NA( I )«NAa * - • ■ • . 

' 59 GO TO 90 

— 55 103 NI2)“1 ' - 

56 N 'A ( 2 ) * N A 8 





SJ GO TO 90 

— 58— -lOH M3>al 

59 NM3J«NAB 

r 60 - -- GO TO 90 

' 61 lOS N(9)>1 

62 - hrf(H)>:Ni<8 - - 

61 ' GO TO 70 

,Q^ N<5J«I 

65 N6 i 5 ] oi^yjS 

6ft - ■ GO TO 90 - - • 

67 107 N(61«l 

68 - 6] ■Nft'S 

69 GO TO 90 

- 7 0 j03 M7>«i . 

71 NAI7)*Na8 

-72 - 60 TO 90 

73 1_09 N(a)«I 

7*t ■ h>M 6}KM.t6 » • 

75 60 TO 90 

74 2,00 CONTtWUe 

77 IF tT.GT.O.COl) 60 TO 21S 

- 76 1F()((‘)95).LT*0*5)G0T0210 

79 IF 7 J •NE.'t69) GO 7f) 210 

— SO ' NA ( 7 > = fJ,vTU 

81 N 171=575 

62 210 nRlTE (4.2051 • ( >^ ( 1 1 » 1 » 1 , 8 ) 

83 205 FORMATC///» SECS ♦ ,8 ( 2X . a6 1 / • 599*.falS/‘ 

8*( 

85 215 DO 220 1=1.9 

66 • ■ - - K"N ( I ) 

87 A(I.S1«X(K) 

■88 - 220 A t I ,6 1 « ( A < 1 , 1 >*A « 1 . 2 ) *A ( 1 , 3 A( I . 8 ) *A ( I .5 1 1 75.0 

89 X (570)«A t V ,6 1 

90 - • • — IF (M7).E0.S73 .ANO. T.LT.Hl.l A(7,6)"0. 

91 C ftRITE(6.300)T,(6t I ,S> .l»l .a) 

92 - - IF ( (T-PTl.LT.i.i .06. A'lOO < T , AflS ( X ( 5 9 9 1 > ) . 6T . 1 . ) 60 TO 310 

93 l-P»T 

9q • PT«LP ^ . 

95 ,.»IT£(6,30C)'PT,(A(l,6j,I«l,6) 

9ft — 300 format <f7.1.8F6.3) 

97 IF ( X (5991.GT.C. 1 60 TO 310 

98 - dimension NXP ( 2qJ .RNXP ( 25) 

99 Data NXP/M5.M8,S6S,q9i .57| , I2q. 123| 139* 138, 1S0» 121 . 

100 6 201 .203. 227*221 .233. 299. 260.91 .576, 130. 4<(*2HZ(20b/ 

101 DO 301 1=1,29 

1C2' - J = tjXP< I 1 

103 301 RN!(Ptn«X(J> 

10ft ARITE (6.3C5) RKXP 

lOS C ARI7E (t>.3C5) < X ( I 1 , l = I , 3 2 1 

10*.. 30S format (7X,ftEl6.6» 

107 310 OQ 320 J« 1 .9 

ioa DO 320 !«1.9 

109 320 A ( 1 . J )>A U .W*!! 

HO RETURN 

111 END 
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6 
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■ " 9 


•c«* 


10 

— 11- 
12 

' 

14 
- • 15 
16 

17 

13 

, t Q . 

C«* 

c»* 

20 

— 21 
22 
• - 23 
24 

c«« 

c*# 

Zb 


26 


• 27 

— 

26 

c»* 

- 29 

■ 

30 

c** 


33 

34 

35 - 

36 
37- 
33 

39- - 
40 

•41 -- 

42 

43 -- 

44 

45 
4A 
47 
43 

-4 9 — 
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BLOCK data 
COHHON/STATE/A C 100) 

COHMON/5TATE/B C50) 

COHKON/ ST ATE/C t bO ) 

• C0MM0H/STATE/0(60) , 

COMMON/stATE/E t bO) 

C0MM0N/STATE/F(201 

COMMON/STATE/C(260) 

state 

DATA A/B9t9,2l5«8i3S*7»2]4>0>7*7i I0«6«27»4»6»5>0i iO«i 1- 10 

1 0«*3«25.0*il3>6.52>3,&2«4i4«lii26.0>2C5>t2.5, 1 1 <• Zq 

2 35B.2.253<5«36 4 H iZOtStSl •ZtBS* 4>3* 1 i 2B4*2>3>2’> 120*7i 21> 30 

3 ) 7« I «43«7>5*0* <0«0> l7«0> •50*0*/ Ol-lOO 

FLOW 

DATA a/39»0..204.,224.,2CS..2ll.*212.,2l4.,21S*,231.tO»»0*iO«/ 101-lSO 

COHP 

Data C/4*0»tl«2il>7»S.3«.2S>2«0«i2*>2«0**«21»<2«>B»«3»3«96i3M‘t»i6 lSl-170 

1 i3*Q. |2 , ,3996 i5.3 • .9058. 9.S?» 1 «50S , 6.0H7 » 171-18C 

2 .2224,2,5)7.5*0, ,«3» I 2»CI«7 Idl-200 

PRES 

DATA 0/40*0. ,0.001 .7*n,, 2*90./ 2Cl-2bo 

RES . 

DATA E/3*. 007506, .004. .01502 , . 052S5 , , 0 1 502 2 , 2* 0 . , *01200 » 

1 •040Q,O...O340Q00,.034n,.03Qo3,4.5oS,.45S5,.07b08..0750fl,.02l02. 

- 2.00733 , .007503 , .01502 ,.t,. 03373, 3. 431. .3754, ,00 4 30 2, 15.39, 2, 35 - 

3 .34*5345.. 2252*. 5255,«3003.. 01502. .4&045.2* 744. *6494, 0..0.* 

4 9*0. ,5. 15/ 

INRT 

DATA F/.G00750S,.002.7*0.,.004,.004,0.,.OQ4,,00626,U*0,/ 301-320 

«1SC 

- OATA 6/48*0. .5000. »5oAo. .30* ,0. ,30. »0. ,05. , 1 5. ,400* .61 .6, 2*0» » 321-300 

1 0., 90. 5.0. ,43.5, 5. 194. 30. .30. , 100. . 188. ,40., 381-390 

-2 3*Q..5o.*Bo.,509..28«.562..0.0.37B., 391-400 

3 SO. , 1 50. » 1 R*0» . . 401-420 

- • 34*0, .32*0. , 4*0. .83. , 2.55 .3.60 .9999, ,0, . 421-495 

• 0. ,0*0 .0. ,2» ,0. ,7*0. ,«7. .2*0. , 496-510 

43..1C.,t0..16.i6.>|o..2*0..ta..t4.. 511-520 

5 I 4, ,2. ,-7 . .-1,4 . .0, ,0. ,«i 4. . 1 3*0 . , 521-540 

- • p. ,C. . ICG. .60S. .430* ,0. . 1 3*-l • .401 . , 541-560 

6 7 2 . » . 0 9 , 6 » 7 , 5 *0 . , 8 . 3 , 0. » 561-57C 

-7 *0550.0*,— 1*5. C..90*. .48. .633, .19, .36. 46.. 571—580 

8 10. <-*015.33., 14 -C* 0, 00,10*. 0./ 581-600 

— ' CO«HOM/X I 00/NJ 9 1 ,*.61 0 ( , IN IT 

DATA INIT/C/ 

• DATA NA/' hR‘,» CO*,* SV* , * V02DOT* , 

* ■ ' SVST,.* OlA^*.; TILT*,' LE6V'/ 

data n/561 ,563.562,570.567 ,568,575,600,570/ 

C0«fl6f./RINTR/R1M 10) ,R0U7( lOI 

- OATA R1N121/12.6/ - 

End 


25|- 
261- 
2; !• 
281- 
291- 


'260 

27o 

260 

29o 

300 


SPRT.T GE# 
Da6-603H32»SE 



FOR RC'MOI 

Bet.— RC** 

FOB RCS(O) 

rEI. fiCS 

FOR RC610J' 

Ret. Rt‘ ■ ^ 

for rc 7 t 0 ) 

REI, RC7 

for «c6<0) 

REU RC e ■- 

FOR RC910J 

REL RC9 - - - 

FOR RCtOtO) 

REI. RCIQ 

FOR RC t 1 (0> 

REI. RC I 1 — — 

FOR RCIRIO) 

REU RC 1 *» ■ • — • 

FOR RCiStOJ 

set. • RC 1 S 

FOR RC16tOJ 

R e L R Cl 6 

FOR RC17(03 

FOR RC I 9 {Cl 

Ret. RCl 9 — 

FOR RC20C0J 

REU RC20 -- - - 

FOR RC2U01 

REL RCZI ■“ 

FOR RCFltOl 

REL RCFl • 

FOR RCF2(0) 

REL RCF2 — 

FOR RCF3IOJ 

REL RCF3 : 

FOR SSVCnUO) 

REL • SSv6<T • ^ : 

ELT S£A0*r (01 . • ' ' 

ELT C02DAT{C) "T'- ■' — " 

■ELT GhOOAT (0) • ' ' " ' ‘ 

FOR delay (0) ; = — - ' — 

REL RELAY 

FOR FC’JSatO) — 

REL FCMSw 

FOR S*IN ( 0) 

REL S-MN 

ELT abAIO) — - 

REL r. n A 

FOR ■ RC3 (0) 

REL RCl 

FOR XJOlO) 

REL XlO 

FOR BLKOaT(O) 

REL BlKOAT 

REL- EAEC — 

ELT EAECtO) 




FOR CVS<0» 

- R£L- CVS — - 

FOR RC13(01 . 

REL RC13 - • 

FOR RCI3/SIHP(Q] 

FOR GMOOIf'l/SAtOJ — 

For RCJ2/SA101 

FOR SS02"/Sa (0>-' 

ELT RUNGlUl 

FOR T£RG/SA«0) 

FOB CVS/SAim 

FOR C0NTRL/SA(0) - 

EUT RufJEX(O) 

ELT RyMOJ 

FOR AlGO C C J 

REU algo 

FOR SSOZAiO) T 

BEL SS02* •' 

ELT - LtSTlOJ. 

FOR CONTRL/OuiOj 

FOR Cv3/OM0»-" 

FOR GrtOOTM/OLCO) 

FOR TERG(OJ 

REL TERG 

FOR COmT«l<OT 1 

REL COnTRu 

REL RC12 

FOR RC12«0J 

FOR GROOInIOT 

REL GKOUl.M 


- A8S - RUN • - 


esRKPT-PRlNT.S - 


iSRfelNAL ^AgE m 
OF POOR QUALITY 
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